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Abstract 

• , 

ISAJET 7.69 is now available, replacing version 7.64 and several intermediate versions. 
The complete set of 1-loop radiative corrections to sparticle masses has been included 
in the ISAJET SUSY spectrum calculation, according to the formulae given by Pierce et 
ai, Nucl. Phys. B491, 3 (1997). Previous versions of ISAJET had included the logarithmic 
1-loop corrections by freezing out soft mass parameters at a scale equal to their mass. The 
current calculation includes all finite corrections as well. Crucial contributions to the encod- 
ing of these expressions were made by Tadas Krupovnickas. We have also adjusted slightly 
the input DR gauge couplings at M z to coincide with measured central values. Logarith- 
mic threshold corrections to gauge and Yukawa couplings are accounted for as in previous 
ISAJET versions by changing the beta functions in the RGEs as various soft mass thresholds 
are passed. We have also enlarged the set of SUSY RGEs to include running of the Higgs 
field vevs. The net effect of these changes is to modify various sparticle masses by typically 
1-5% from the predictions of ISAJET 7.64. 

The size of the small mass splitting between Xi an d Xi 1S important for AMSB phe- 
nomenology. Rather than extracting this splitting from the general result, for the AMSB 
model we calculate it from the tree and the 1-loop vector boson loop graphs [taken from 
Cheng, Dobrescu and Matchev, Nucl. Phys. B543, 47, (1999)] and use it to set the xt mass. 
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We have modified the RGE solution to use the previous SUSY masses to compute \i 
before computing the new SUSY masses and from them the loop corrections to /i. This 
seems to improve the convergence, e.g., in the focus point, or hyperbolic branch region, of 
the mSUGRA model. The RGE solution still does not converge properly in isolated regions. 

We have also corrected a sign mistake in the A-term boundary conditions for the anomaly- 
mediated SUSY breaking model. 

The new version can be obtained either by anonymous FTP from 

\protect\vrule widthOpt\protect\href {f tp : //f tp .phy .bnl . gov/pub/isajet}{f tp : //f tp .phy .bnl 
or by HTTP from 

\protect\vrule widthOpt\protect\href {http : //www. phy .bnl .gov\string~isajet}{http : //www.ph 

or via AFS from 

/af s/cern . ch/user/p/paige/public/isaj et 

All of these contain the same card image PAM file isajet.car and Unix Makefile to build 
IS A JET from isajet.car on all supported Unix platforms. 
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1 Introduction 



ISAJET is a Monte Carlo program which simulates pp, pp and e + e~ interactions at high 
energies. ISAJET is based on perturbative QCD plus phenomenological models for parton 
and beam jet fragmentation. Events are generated in four distinct steps: 

• A primary hard scattering is generated according to the appropriate QCD cross section. 

• QCD radiative corrections are added for both the initial and the final state. 

• Partons are fragmented into hadrons independently, and particles with lifetimes less 
than about 1CT 12 seconds are decayed. 

• Beam jets are added assuming that these are identical to a minimum bias event at the 
remaining energy. 

ISAJET incorporates ISASUSY, which evaluates branching ratios for the minimal su- 
persymmetric extension of the standard model. H. Baer and X. Tata are coauthors of this 
package, and they have done the original calculations with various collaborators. See the 
ISASUSY documentation in the patch Section 12. 

ISAJET is supported for ANSI Fortran and for Cray, DEC Ultrix, DEC VMS, HP/9000 
7xx, IBM VM/CMS 370 and 30xx, IBM AIX RS/6000, Linux, Silicon Graphics 4D, and Sun 
computers. The CDC 7600 and ETA 10 versions are obsolete and are no longer supported. 
It is written mainly in ANSI standard FORTRAN 77, but it does contain some extensions 
except in the ANSI version. The code is maintained with a combination of RCS, the Revision 
Control System, and the Patchy code management system, which is part of the CERN 
Library. The original sources are kept on physgi01.phy.bnl.gov in "isajet/ isalibrary/RCS; 
decks revised in release n.nn are kept in ~isajet/isalibrary/nnn. ISAJET is supplied to 
BNL, CERN, Fermilab, and SLAC; it is also available from 

\protect\vrule widthOpt\protect\href {http : //www . phy . bnl . gov\string~isajet}{http : //www.ph 
\protect\vrule widthOpt\protect\href {f tp : //ftp .phy .bnl . gov/pub/isajet}{f tp : //ftp .phy .bnl 

or by request from the authors. 

Patch ISAPLT contains the skeleton of an HBOOK histogramming job, a trivial calorime- 
ter simulation, and a jet-finding algorithm. (The default is HBOOK4; HBOOK3 can be 
selected with a Patchy switch.) These are provided for convenience only and are not sup- 
ported. 
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2 Physics 



ISAJET is a Monte Carlo program which simulates pp, pp and e + e~ interactions at high 
energy. The program incorporates perturbative QCD cross sections, initial state and final 
state QCD radiative corrections in the leading log approximation, independent fragmentation 
of quarks and gluons into hadrons, and a phenomenological model tuned to minimum bias 
and hard scattering data for the beam jets. 

2.1 Hard Scattering 

The first step in simulating an event is to generate a primary hard scattering according to 
some QCD cross section. This has the general form 

a = a F{x 1 ,Q 2 )F{x 2 ,Q 2 ) 

where <r is a cross section calculated in QCD perturbation theory, F(x, Q 2 ) is a struc- 
ture function incorporating QCD scaling violations, x\ and x 2 are the usual parton model 
momentum fractions, and Q 2 is an appropriate momentum transfer scale. 

For each of the processes included in ISAJET, the basic cross section a is a two-body 
one, and the user can set limits on the kinematic variables and type for each of the two 
primary jets. For DRELLYAN and WPAIR events the full matrix element for the decay of 
the W's into leptons or quarks is also included. 

The following processes are available: 

2.1.1 Minbias 

No hard scattering at all, so that the event consists only of beam jets. Note that at high 
energy the jet cross sections become large. To represent the total cross section it is better 
to use a sample of TWOJET events with the lower limit on pt chosen to give a cross section 
equal to the inelastic cross section or to use a mixture of MINBIAS and TWOJET events. 

2.1.2 Twojet 

All order a 2 QCD processes, which give rise in lowest order to two high-p t jets. Included 
are, e.g. 

9 + 9 -»• 9 + 9 

g + q -> g + q 
9 + 9 -> q + q 

Masses are neglected for c and lighter quarks but are taken into account for b and t quarks. 
The Q 2 scale is taken to be 

Q 2 = 2stu/(s 2 + t 2 + u 2 ) 

The default parton distributions are those of the CTEQ Collaboration, fit CTEQ5L, using 
lowest order QCD evolution. A few older sets of parton distributions are included. There 
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is also an interface to the CERN PDFLIB compilation of parton distributions. Note that 
structure functions for heavy quarks are included, so that processes like 

g + t^g + t 

can be generated. The Duke-Owens parton distributions do not contain b or t quarks. 

Since the t is so heavy, it decays before it can hadronize, so instead of t hadrons a t 
quark appears in the particle list. It is decayed using the V — A matrix element including 
the W propagator with a nonzero width, so the same decays should be used for writ < mw 
and m t > mw] the W should not be listed as part of the decay mode. The partons are then 
evolved and fragmented as usual; see below. The real or virtual W and the final partons 
from the decay, including any radiated gluons, are listed in the particle table, followed by 
their fragmentation products. Note that for semileptonic decays the leptons appear twice: 
the lepton parton decays into a single particle of the same type but in general somewhat 
different momentum. In all cases only particles with IDCAY = should be included in the 
final state. 

A fourth generation x, y is also allowed. Fourth generation quarks are produced only 
by gluon fusion. Decay modes are not included in the decay table; for a sequential fourth 
generation they would be very similar to the t decays. In decays involving quarks, it is 
essential that the quarks appear last. 

2.1.3 Drellyan 

Production of a W in the standard model, including a virtual 7, a W + , a W~, or a Z°, and 
its decay into quarks or leptons. If the transverse momentum QTW of the W is fixed equal 
to zero then the process simulated is 

q + q —> W — > q + q 
-> £ + £ 

Thus the W has zero transverse momentum until initial state QCD corrections are taken 
into account. If non-zero limits on the transverse momentum q t for the W are set, then 
instead the processes 

q + q -> W + g 
g + q -> W + q 

are simulated, including the full matrix element for the W decay. These are the dominant 
processes at high q t) but they are of course singular at q t = 0. A cutoff of the l/g t 2 singularity 
is made by the replacement 

1/^1/^+4 ql = (.2 GeV)M 

This cutoff is chosen to reproduce approximately the qt dependence calculated by the sum- 
mation of soft gluons and to give about the right integrated cross section. Thus this option 
can be used for low as well as high transverse momenta. 
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The scale for QCD evolution is taken to be proportional to the mass for lowest order 
Drell-Yan and to the transverse momentum for high-p t Drell-Yan. The constant is adjusted 
to get reasonable agreement with the W + n jet cross sections calculated from the full QCD 
matrix elements by F.A. Berends, et al., Phys. Lett. B224, 237 (1989). 

For the processes g + b — > W + t and g + t — > Z + t, cross sections with a non-zero top 
mass are used for the production and the W/Z decay. These were calculated using FORM 
1.1 by J. Vermaseren. The process g + t ^ W + b is not included. Both g + b — > W~ + t and 
g + t — > W~ + 6 of course give the same W~ + t + b final state after QCD evolution. While 
the latter process is needed to describe the m t = 0(!) mass singularity for q t ^> m t , it has a 
pole in the physical region at low qt from on-shell t W + b decays. There is no obvious 
way to avoid this without introducing an arbitrary cutoff. Hence, selecting only W + b will 
produce a zero cross section. The Q 2 scale for the parton distributions in these processes 
is replaced by Q 2 + m 2 ; this seems physically sensible and prevents the cross sections from 
vanishing at small q t . 

2.1.4 Photon 

Single and double photon production through the lowest order QCD processes 

g + q -> 7 + <? 
q + q -> 7 + 3 
q + q -> 7 + 7 

Higher order corrections are not included. But 7's, W's, and Z's are radiated from final 
state quarks in all processes, allowing study of the bremsstrahlung contributions. 

2.1.5 Wpair 

Production of pairs of W bosons in the standard model through quark-antiquark annihilation, 

q + q -> W + + W~ 
-> Z° + Z° 

-»• W^ + + 7,H/-+7 
- ^° + 7 

The full matrix element for the W decays, calculated in the narrow resonance approximation, 
is included. However, the higher order processes, e.g. 

q + q-^q + q + W + + W~ 

are ignored, although they in fact dominate at high enough mass. Specific decay modes can 
be selected using the WMODEi keywords. 
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2.1.6 Higgs 

Production and decay of the standard model Higgs boson. The production processes are 

g + g — > H (through a quark loop) 

q + q — > H (with t + i dominant) 

W + + W~ — > if (with longitudinally polarized W) 

Z° + Z° -> if (with longitudinally polarized Z) 

If the (Standard Model) Higgs is lighter than 2M W) then it will decay into pairs of fermions 
with branching ratios proportional to mj. If it is heavier than 2Mw, then it will decay 
primarily into W + W~ and Z°Z° pairs with widths given approximately by 



r(H -> z°z°) 

Numerically these give approximately 

= 0.5 TeV 



G F Mfj 
8tt^2 
G F Mfj 



M H ^ 



3 



ITeV 



The width proportional to arises from decays into longitudinal gauge bosons, which like 
Higgs bosons have couplings proportional to mass. 

Since a heavy Higgs is wide, the narrow resonance approximation is not valid. To obtain 
a cross section with good high energy behavior, it is necessary to include a complete gauge- 
invariant set of graphs for the processes 



w + w~ 


-> w + w 


w + w~ 




z°z° 




z°z° 





with longitudinally polarized H^~, and Z° bosons in the initial state. This set of 

graphs and the corresponding angular distributions for the W + , W~, and Z° decays have 
been calculated in the effective W approximation and included in HIGGS. The W structure 
functions are obtained by integrating the EHLQ parameterization of the quark ones term by 
term. The Cabibbo-allowed branchings 

q -> W + + q' 
q -> W~ + q' 
q -> Z° + q 

are generated by backwards evolution, and the standard QCD evolution is performed. This 
correctly describes the W collinear singularity and so contains the same physics as the 
effective W approximation. 
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If the Higgs is lighter than 2M W , then its decay to 77 through W and t loops may be 
important. This is also included in the HIGGS process and may be selected by choosing GM 
as the jet type for the decay. 

If the Higgs has Mz < Mh < 2Mz, then decays into one real and one virtual Z° are 
generated if the ZO ZO decay mode is selected, using the calculation of Keung and Marciano, 
Phys. Rev. D30, 248 (1984). Since the calculation assumes that one Z° is exactly on shell, 
it is not reliable within of order the Z° width of M H = 2M Z ; Higgs and and Z°Z° masses 
in this region should be avoided. The analogous Higgs decays into one real and one virtual 
charged W are not included. 

Note that while HIGGS contains the dominant graphs for Higgs production and graphs 
for W pair production related by gauge invariance, it does not contain the processes 

q + q -> W + W~ 
q + q -> Z°Z° 

which give primarily transverse gauge bosons. These must be generated with WPAIR. 

If the MSSMi or SUGRA keywords are used with HIGGS, then one of the three MSSM neutral 
Higgs is generated instead using gluon-gluon and quark-antiquark fusion with the appropriate 
SUSY couplings. Since heavy CP even SUSY Higgs are weakly coupled to W pairs and CP 
odd ones are completely decoupled, WW fusion and WW — > WW scattering are not included 
in the SUSY case. (WW — > WW can be generated using the Standard Model process with 
a light Higgs mass, say 100 GeV.) The MSSM Higgs decays into both Standard Model and 
SUSY modes as calculated by ISASUSY are included. For more discussion see the SUSY 
subsection below and the writeup for ISASUSY. The user must select which Higgs to generate 
using HTYPE; see Section 6 below. If a mass range is not specified, then the range mass 
M H ± hY H is used by default. (This cannot be done for the Standard Model Higgs because 
it is so wide for large masses.) Decay modes may be selected in the usual way. 

2.1.7 WHiggs 

Generates associated production of gauge and Higgs bosons, i.e., 

q + q^H + W,H + Z, 

in the narrow resonance approximation. The desired subprocesses can be selected with 
JETTYPEi, and specific decay modes of the W and/or Z can be selected using the WMODEi 
keywords. Standard Model couplings are assumed unless SUSY parameters are specified, in 
which case the SUSY couplings are used. 

2.1.8 SUSY 

Generates pairs of supersymmetric particles from gluon-quark or quark-antiquark fusion. If 
the MSSMi or SUGRA parameters defined in Section 6 below are not specified, then only 
gluinos and squarks are generated: 

9 + 9 -> 9 + 9 
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g + g -> g + g 
g + q -> # + g 

# + # -> <? + g~ 
g + g — ► g + g 
g + g — > g + g 

Left and right squarks are distinguished but assumed to be degenerate. Masses can be 
specified using the GAUGINO, SQUARK, and SLEPTON parameters described in Section 6. No 
decay modes are specified, since these depend strongly on the masses. The user can either 
add new modes to the decay table (see Section 9) or use the FORCE or F0RCE1 commands 
(see Section 6). 

If MSSMA, MSSMB, and MSSMC are specified, then the ISASUSY package is used to calculate 
the masses and decay modes in the minimal supersymmetric extension of the standard model 
(MSSM), assuming SUSY grand unification constraints in the neutralino and chargino mass 
matrix but allowing some additional flexibility in the masses. The scalar particle soft masses 
are input via MSSMi, so that the physical masses will be somewhat different due to JJ-term 
contributions and mixings for 3rd generation sparticles. t\ and ti production and decays are 
now included. The lightest SUSY particle is assumed to be the lightest neutralino Z\. If the 
MSSMi parameters are specified, then the following additional processes are included using 
the MSSM couplings for the production cross sections: 

g + q -> Zi + q, W { + q 

q + q -> Zi + g, W t + g 

q + q -> Wi + Zj 

q + q -> W+ + W- 

q + q -> Zi + Zj 

q + q -> £ + + £~, v + v 

Processes can be selected using the optional parameters described in Section 6 below. 

Beginning with Version 7.42, matrix elements are taken into account in the event gener- 
ator as well as in the calculation of decay widths for MSSM three-body decays of the form 
A — > Bff, where A and B are gluinos, charginos, or neutralinos. This is implemented 
by having ISASUSY save the poles and their couplings when calculating the decay width 
and then using these to reconstruct the matrix element. Other three-body decays may be 
included in the future. Decays selected with FORCE use the appropriate matrix elements. 

An optional keyword MSSMD can be used to specify the second generation masses, which 
otherwise are assumed degenerate with the first generation. An optional keyword MSSME can 
be used to specify values of the U(l) and 577(2) gaugino masses at the weak scale rather 
than using the default grand unification values. The chargino and neutralino masses and 
mixings are then computed using these values. 
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2.1.9 SUSY Models 



The 24 MSSMi parameters describe the MSSM at the weak scale with the additional assump- 
tions of exact flavor and CP conservation; the general MSSM has 105 parameters. These 
weak-scale SUSY-breaking parameters presumably arise from spontaneous SUSY breaking 
in a hidden sector that is communicated to the MSSM at some scale M 3> M z . There are a 
number of plausible models in which this symmetry breaking is simple, so that the MSSM at 
the high scale involves only a small number of parameters. These are then related to those 
at the weak scale by the renormalization group equations (RGE's). 

Isajet therefore implements in subroutine SUGRA the complete 2-loop RGE's for the 
gauge couplings, Yukawa couplings, and soft breaking terms. Contributions from right- 
handed neutrinos are optionally included. The RGE's are solved iteratively, running from the 
weak scale to the high scale M and back using Runge-Kutta integration. After each iteration 
the SUSY masses are recalculated, and the renormalization group improved 1-loop corrected 
Higgs potential is calculated and minimized. These results are used to modify the RGE 
/9-functions appropriately as each threshold is crossed during the next iteration. Beginning 
with Version 7.65, the constant parts as well as the logarithms from these thresholds are 
included using the results of Pierce et al, Nucl. Phys. B491, 3 (1997). The whole process is 
repeated, increasing the number of Runge-Kutta steps by a factor of 1.2 for each iteration, 
until all the RGE variables except fi and B differ by less than 0.3%. Since /i and B vary 
rapidly near the weak scale, they are only required to differ by less than 5%. The requirement 
of good electroweak symmetry breaking, /i 2 > 0, is only imposed after the iterative solution 
has converged. 

A number of different models for SUSY breaking at the high scale are included in ISAJET. 
The SUGRA parameters must be specified for the minimal supergravity framework. This 
assumes that the gauge couplings unify at the GUT scale, M ~ 10 16 GeV, defined by a.2 = 
a\. SUSY breaking occurs at that scale with universal soft breaking terms produced by 
gravitational interactions with a hidden sector. The parameters of the model are 

• m : the common scalar mass at the GUT scale; 

• mi/ 2 : the common gaugino mass at the GUT scale; 

• A : the common soft trilinear SUSY breaking parameter at the GUT scale; 

• tan/5: the ratio of Higgs vacuum expectation values at the electroweak scale; 

• sgn^u = ±1: the sign of the Higgsino mass term. 

An attractive feature of this model is that the Higgs are unified with the other scalars at the 
GUT scale but m 2 Hu is driven negative by the large top Yukawa coupling f t . Isajet imposes 
this radiative symmetry breaking for the SUGRA model but not other possible constraints 
such as b-r unification or limits on proton decay. 

The SUGRA model with exact compling constant unification produces too large a value 
of a s at the weak scale. The default is to use the experimental value, assuming that threshold 
effects at the GUT scale produce this. Exact unification can also be imposed. 
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The assumption of universality at the GUT scale is rather restrictive and may not be 
valid. A variety of non-universal SUGRA (NUSUGRA) models can be generated using the 
NUSUG1, . . . , NUSUG5 keywords. These might be used to study how well one could test the 
minimal SUGRA model. The keyword SSBCSC can be used to specify an alternative scale 
(i.e., not the coupling constant unification scale) for the RGE boundary conditions. 

An alternative to the SUGRA model is the Gauge Mediated SUSY Breaking (GMSB) 
model of Dine, Nelson, and collaborators. In this model SUSY breaking is communicated 
through gauge interactions with messenger fields at a scale M m small compared to the Planck 
scale and are proportional to gauge couplings times A m . Since M m is small and the masses 
at it are the same for each generation, there are no flavor changing neutral currents. The 
messenger fields should form complete SU(5) representations to preserve the unification of 
the coupling constants. The parameters of the GMSB model, which are specified by the 
GMSB keyword, are 

• A m = F m /M m : the scale of SUSY breaking, typically 10-100 TeV; 

• M m > A m : the messenger mass scale; 

• N 5 : the equivalent number of 5 + 5 messenger fields. 

• tan/5: the ratio of Higgs vacuum expectation values at the electroweak scale; 

• sgn/i = ±1: the sign of the Higgsino mass term; 

• Cgrav > 1: the ratio of the gravitino mass to the value it would have had if the only 
SUSY breaking scale were F m . 

In GMSB models the lightest SUSY particle is always the nearly massless gravitino G. The 
parameter C grav scales the gravitino mass and hence the lifetime of the next lightest SUSY 
particle to decay into it. The NOGRAV keyword can be used to turn off gravitino decays. 

A variety of non-minimal GMSB models can be generated using additional parameters 
set with the GMSB2 keyword. These additional parameters are 

• ft, an extra factor multiplying the gaugino masses at the messenger scale. (Models 
with multiple spurions generally have ft < 1.) 

• 5Mfj d1 SMft u , Higgs mass-squared shifts relative to the minimal model at the messenger 
scale. (These might be expected in models which generate \x realistically.) 

• Dy[M\ a U(l)y messenger scale mass-squared term (D-term) proportional to the 
hyper charge Y . 

• JV 5l , N 52 , and N 5s , independent numbers of gauge group messengers. They can be 
non-integer in general. 

For discussions of these additional parameters, see S. Dimopoulos, S. Thomas, and J.D. 
Wells, hep-ph/9609434, Nucl. Phys. B488, 39 (1997), and S.P Martin, hep-ph/9608224, 
Phys. Rev. D55, 3177 (1997). 
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Gravitino decays can be included in the general MSSM framework by specifying a grav- 
itino mass with MGVTNO. The default is that such decays do not occur. 

Another alternative SUSY model choice allowed is anomaly-mediated SUSY breaking, 
developed by Randall and Sundrum. In this model, it is assumed that SUSY breaking takes 
place in other dimensions, and SUSY breaking is communicated to the visible sector via the 
superconformal anomaly. In this model, the lightest SUSY particle is usually the neutralino 
which is nearly pure wino-like. The chargino is nearly mass degenerate with the lightest 
neutralino. It can be very long lived, or decay into a very soft pion plus missing energy. 
The model incorporated in ISAJET, based on work by Ghergetta, Giudice and Wells (hep- 
ph/9904378), and by Feng and Moroi (hep-ph/9907319) adds a universal contribution Wq 
to all scalar masses to avoid problems with tachyonic scalars. The parameters of the model, 
which can be set via the AMSB keyword, are 

• m : Common scalar mass; 

• m 3 / 2 : Gravitino mass, typically lOTeV since rrii = (Pi/ gi)m 3 / 2 - 

• tan/?: Usual ratio of vev's at weak scale; 

• sgn/i: Usual sign of fi, ±1. 

Care should be taken with the chargino decay, since it may have macroscopic decay lengths, 
or even decay outside the detector. 

Since neutrinos seem to have mass, the effect of a massive right-handed neutrino has been 
included in ISAJET, when calculating the sparticle mass spectrum. If the keyword SUGRHN 
is used, then the user must input the 3rd generation neutrino mass (at scale Mz) in units 
of GeV, and the intermediate scale right handed neutrino Majorana mass Mjv, also in GeV. 
In addition, one must specify the soft SUSY-breaking masses A n and mo R valid at the GUT 
scale. Then the neutrino Yukawa coupling is computed in the simple see-saw model, and 
renormalization group evolution includes these effects between Mqut and M^. Finally, to 
facilitate modeling of SO (10) SUSY-GUT models, loop corrections to 3rd generation fermion 
masses have been included in the ISAJET SUSY models. 

The ISASUSY program can also be used independently of the rest of ISAJET, either 
to produce a listing of decays or in conjunction with another event generator. Its physics 
assumptions are described in more detail in Section 12. ISASUSY accepts soft SUSY breaking 
parameters at the weak scale and calculates the masses and decay modes from them. The 
ISASUGRA program can also be used independently to solve the renormalization group 
equations with SUGRA, NUSUGRA, GMSB, or AMSB boundary conditions and then to 
call ISASUSY to calculate the decay modes. The main programs described in Section 5.4 
prompt for interactive input and print the results to a file. 

Generally the MSSM, SUGRA, or GMSB option should be used to study supersymmetry 
signatures; the SUGRA or GMSB parameter space is clearly more manageable. The more 
general option may be useful to study alternative SUSY models. It can also be used, e.g., to 
generate pointlike color-3 leptoquarks in technicolor models by selecting squark production 
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and setting the gluino mass to be very large. The MSSM or SUGRA option may also be 
used with top pair production to simulate top decays to SUSY particles. 



2.1.10 e+e" 

An e + e~ event generator is also included in ISAJET. The Standard Model processes included 
are e + e~ annihilation through 7 and Z to quarks and leptons, and production of W + W~ 
and Z°Z° pairs. In contrast to WPAIR and HIGGS for the hadronic processes, the produced 
W's and Z's are treated as particles, so their spins are not properly taken into account in 
their decays. (Because the Ws and Z's are treated as particles, their decay modes can be 
selected using FORCE or F0RCE1, not WMODEi. See Section [6] below.) Other Standard Model 
processes, including e + e~ — > e + e~ (t-channel graph) and e + e~ — > 77, are not included. Once 
the primary reaction has been generated, QCD radiation and hadronization are done as for 
hadronic processes. 

The e + e~ generator can be run assuming no initial state radiation (the default), or 
an initial state electron structure function can be used for bremsstrahlung or the combina- 
tion bremsstrahlung/beamstrahlung effect. Bremsstrahlung is implemented using the Fadin- 
Kuraev e~ distribution function, and can be turned on using the EEBREM command while 
stipulating the minimal and maximal subprocess energy. Beamstrahlung is implemented by 
invoking the EEBEAM keyword. In this case, in addition the beamstrahlung parameter T and 
longitudinal beam size a z (in mm) must be given. The definition for T in terms of other beam 
parameters can be found in the article Phys. Rev. D49, 3209 (1994) by Chen, Barklow and 
Peskin. The bremsstrahlung structure function is then convoluted with the beamstrahlung 
distribution (as calculated by P. Chen) and a spline fit is created. Since the cross section 
can contain large spikes, event generation can be slow if a huge range of subprocess energy 
is selected for light particles; in these scenarios, NTRIES must be increased well beyond the 
default value. 

e + e~ annihilation to SUSY particles is included as well with complete lowest order dia- 
grams, and cascade decays. The processes include 

e + e~ — > qq 

e + e~ -> it 

e + e~ -> WiWj 

e + e~ — > ZiZj 

e + e- -> H L + Z ,H H + Z ,H A + H L ,H A + H° H ,H + + H- 

Note that SUSY Higgs production via WW and ZZ fusion, which can dominate Higgs 
production processes at y/s > 500 GeV, is not included. Spin correlations are neglected, 
although 3-body sparticle decay matrix elements are included. 

e + e~ cross sections with polarized beams are included for both Standard Model and 
SUSY processes. The keyword EPOL is used to set Pl{g~) and P L (e + ), where 

P L (e) = (n L - n R )/ (n L + n R ) 

so that — 1 < Pl < +1. Thus, setting EPOL to —.9, will yield a 95% right polarized electron 
beam scattering on an unpolarized positron beam. 
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2.1.11 Technicolor 



Production of a technirho of arbitrary mass and width decaying into W Z° or W + W~ pairs. 
The cross section is based on an elastic resonance in the WW cross section with the effective 
W approximation plus a W mixing term taken from EHLQ. Additional technicolor processes 
may be added in the future. 

2.1.12 Extra Dimensions 

The possibility that there might be more than four space-time dimensions at a distance scale 
R much larger than G]!, 2 has recently attracted interest. In these theories, 

C 1 

N 8tcR*M 2 d +s ' 

where 5 is the number of extra dimensions and Mo is the 4 + 5 Planck scale. Gravity 
deviates from the standard theory at a distance R ~ 10 22 / 5 - 19 m, so 5 > 2 is required. If 
M D is of order 1 TeV, then the usual heirarchy problem is solved, although there is then a 
new heirarchy problem of why R is so large. 

In such models the graviton will have many Kaluza-Klein excitations with a mass splitting 
of order 1/R. While any individual mode is suppressed by the four-dimensional Planck mass, 
the large number of modes produces a cross section suppressed only by l/M D . The signature 
is an invisible massive graviton plus a jet, photon, or other Standard Model particle. The 
EXTRADIM process implements this reaction using the cross sections of Giudice, Rattazzi, 
and Wells, hep-ph/9811291. The number 5 of extra dimensions, the mass scale Md, and the 
logical flag UVCUT are specified using the keyword EXTRAD. If UVCUT is TRUE, the cross section 
is cut off above the scale M D ; the model is not valid if the results depend on this flag. 

2.2 Multiparton Hard Scattering 

All the processes listed in Section 2.1 are either 2^2 processes like TWO JET or 2 — > 1 
s-channel resonance processes followed by a 2-body decay like DRELLYAN. The QCD parton 
shower described in Section 2.3 below generates multi-parton final states starting from these, 
but it relies on an approximation which is valid only if the additional partons are collinear 
either with the initial or with the final primary ones. Since the QCD shower uses exact 
non-colliear kinematics, it in fact works pretty well in a larger region of phase space, but it 
is not exact. 

Non-collinear multiparton final states are interesting both in their own right and as back- 
grounds for other signatures. Both the matrix elements and the phase space for multiparton 
processes are complicated; they have been incorporated into ISAJET for the first time in 
Version 7.45. To calculate the matrix elements we have used the MadGraph package by 
Stelzer and Long, Comput. Phys. Commun. 81, 357 (1994), hep-ph/9401258. This auto- 
matically generates the amplitude using HELAS, a formalism by Murayama, Watanabe, and 
Hagiwarak KEK-91-11, that calculates the amplitude for any Feynman diagram in terms of 
spinnors, vertices, and propagators. The MadGraph code has been edited to incorporate 
summations over quark flavors. To do the phase space integration, we have used a simple 
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recursive algorithm to generate n-body phase space. We have included limits on the total 
mass of the final state using the MTOT keyword. Limits on the p T and rapidity of each final 
parton can be set via the PT and Y keyworks, while limits on the mass of any pair of final 
partons can be set via the MIJTOT keyword. These limits are sufficient to shield the infrared 
and collinear singularities and to render the result finite. However, the parton shower pop- 
ulates all regions of phase space, so careful thought is needed to combine the parton-shower 
based and multiparton based results. 

While the multiparton formalism is rather general, it still takes a substantial amount of 
effort to implement any particular process. So far only one process has been implemented. 

2.2.1 Z + 2jets 

The ZJJ process generates a Z boson plus two jets, including the qq — > Zqq, gg — > Zqq, 
qq — > Zgg, qq — > Zqq, and gq — > Zgq processes. The Z is defined to be jet 1; it is treated in 
the narrow resonance approximation and is decayed isotropically. The quarks, antiquarks, 
and gluons are defined to be jets 2 and 3 and are symmetrized in the usual way. 

2.3 QCD Radiative Corrections 

After the primary hard scattering is generated, QCD radiative corrections are added to allow 
the possibility of many jets. This is essential to get the correct event structure, especially at 
high energy. 

Consider the emission of one extra gluon from an initial or a final quark line, 

q(p) -> q(p 1 ) + g(p 2 ) 

From QCD perturbation theory, for small p 2 the cross section is given by the lowest order 
cross section multiplied by a factor 

a = a a s (p 2 )/(2np 2 )P(z) 

where z — Pi/p and P(z) is an Altarelli-Parisi function. The same form holds for the other 
allowed branchings, 

g(p) -> g(pi)+g(p 2 ) 
g(p) -> q{p l ) + q(p 2 ) 

These factors represent the collinear singularities of perturbation theory, and they produce 
the leading log QCD scaling violations for the structure functions and the jet fragmentation 
functions. They also determine the shape of a QCD jet, since the jet M 2 is of order a s p\ 
and hence small. 

The branching approximation consists of keeping just these factors which dominate in the 
collinear limit but using exact, non-collinear kinematics. Thus higher order QCD is reduced 
to a classical cascade process, which is easy to implement in a Monte Carlo program. To 
avoid infrared and collinear singularities, each parton in the cascade is required to have a 
mass (spacelike or timelike) greater than some cutoff t c . The assumption is that all physics 
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at lower scales is incorporated in the nonperturbative model for hadronization. In ISAJET 
the cutoff is taken to be a rather large value, (6 GeV) 2 , because independent fragmentation 
is used for the jet fragmentation; a low cutoff would give too many hadrons from overlapping 
partons. It turns out that the branching approximation not only incorporates the correct 
scaling violations and jet structure but also reproduces the exact three-jet cross section 
within factors of order 2 over all of phase space. 

This approximation was introduced for final state radiation by Fox and Wolfram. The 
QCD cascade is determined by the probability for going from mass t to mass ti emitting 
no resolvable radiation. For a resolution cutoff z c < z < 1 — z c , this is given by a simple 
expression, 

P(t ,t 1 ) = (a s (t )/a s (t 1 )) 2 ^ )/b0 

where 

l{z c ) = [ 1ZC dz P(z), b = (33 - 2n / )/(127r) 

Clearly if P(t ,ti) is the integral probability, then dP/dti is the probability for the first 
radiation to occur at t\. It is straightforward to generate this distribution and then iteratively 
to correct it to get a cutoff at fixed t c rather than at fixed z c . 

For the initial state it is necessary to take account of the spacelike kinematics and of the 
structure functions. Sjostrand has shown how to do this by starting at the hard scattering 
and evolving backwards, forcing the ordering of the spacelike masses t. The probability 
that a given step does not radiate can be derived from the Altarelli-Parisi equations for the 
structure functions. It has a form somewhat similar to P(to, t\) but involving a ratio of the 
structure functions for the new and old partons. It is possible to find a bound for this ratio 
in each case and so to generate a new t and z as for the final state. Then branchings for 
which the ratio is small are rejected in the usual Monte Carlo fashion. This ratio suppresses 
the radiation of very energetic partons. It also forces the branching g — > t + i for a t quark 
if the t structure function vanishes at small momentum transfer. 

At low energies, the branching of an initial heavy quark into a gluon sometimes fails; 
these events are discarded and a warning is printed. 

Since t c is quite large, the radiation of soft gluons is cut off. To compensate for this, 
equal and opposite transverse boosts are made to the jet system and to the beam jets after 
fragmentation with a mean value 

( P 2 t ) = (.lGeV)^ 

The dependence on Q 2 is the same as the cutoff used for DRELLYAN and the coefficient is 
adjusted to fit the p t distribution for the W. 

Radiation of gluons from gluinos and scalar quarks is also included in the same approxi- 
mation, but the production of gluino or scalar quark pairs from gluons is ignored. Very little 
radiation is expected for heavy particles produced near threshold. 

Radiation of photons, W's, and Z's from final state quarks is treated in the same approxi- 
mation as QCD radiation except that the coupling constant is fixed. Initial state electroweak 
radiation is not included; it seems rather unimportant. The W + 's, W~'s and Z's are de- 
cayed into the modes allowed by the WPMODE, WMMODE, and ZOMODE commands respectively. 
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Warning: The branching ratios implied by these commands are not included in the cross 
section because an arbitrary number of W's and Z's can in principle be radiated. 

2.4 Jet Fragmentation: 

Quarks and gluons are fragmented into hadrons using the independent fragmentation ansatz 
of Field and Feynman. For a quark q, a new quark- ant iquark pair q^ is generated with 

u : d : s = .43 : .43 : .14 

A meson qqi is formed carrying a fraction z of the momentum, 

E' + p' z = z(E + Pz ) 

and having a transverse momentum p t with (p t ) = 0.35 GeV. Baryons are included by 
generating a diquark with probability 0.10 instead of a quark; adjacent diquarks are not 
allowed, so no exotic mesons are formed. For light quarks z is generated with the splitting 
function 

f(z) = 1 - a + a(b + 1)(1 - z)\ a = 0.96, b = 3 
while for heavy quarks the Peterson form 

f(z)=x(l-xf/((l-xf + exf 

is used with e = .80/m 2 for c and e = .50/m 2 for q = b,t,y,x. These values of e have 
been determined by fitting PEP, PETRA, and LEP data with ISAJET and should not be 
compared with values from other fits. Hadrons with longitudinal momentum less than zero 
are discarded. The procedure is then iterated for the new quark q\ until all the momentum 
is used. A gluon is fragmented like a randomly selected u, d, or s quark or antiquark. 

In the fragmentation of gluinos and scalar quarks, supersymmetric hadrons are not dis- 
tinguished from partons. This should not matter except possibly for very light masses. The 
Peterson form for f(x) is used with the same value of epsilon as for heavy quarks, e = 0.5/m 2 . 

Independent fragmentation correctly describes the fast hadrons in a jet, but it fails to 
conserve energy or flavor exactly. Energy conservation is imposed after the event is generated 
by boosting the hadrons to the appropriate rest frame, rescaling all of the three-momenta, 
and recalculating the energies. 

2.5 Beam Jets 

There is now experimental evidence that beam jets are different in minimum bias events and 
in hard scattering events. ISAJET therefore uses similar a algorithm but different parameters 
in the two cases. 

The standard models of particle production are based on pulling pairs of particles out 
of the vacuum by the QCD confining field, leading naturally to only short-range rapidity 
correlations and to essentially Poisson multiplicity fluctuations. The minimum bias data 
exhibit KNO scaling and long-range correlations. A natural explanation of this was given 
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by the model of Abramovskii, Kanchelli and Gribov. In their model the basic amplitude 
is a single cut Pomeron with Poisson fluctuations around an average multiplicity (n), but 
unitarity then produces graphs giving K cut Pomerons with multiplicity K(n). 

A simplified version of the AKG model is used in ISAJET. The number of cut Pomerons is 
chosen with a distribution adjusted to fit the data. For a minimum bias event this distribution 
is 

P{K) = (1 + AK 2 ) exp -1.8K 

while for hard scattering 

P(l) -> O.LP(l), P(2) -> 0.2P(2), P(3)^0.5P(3) 

For each side of each event an x for the leading baryon is selected with a distribution varying 
from flat for K = 1 to like that for mesons for large K: 

f(x) = N(K)(1 - xo) c (K), c(K) = l/K + (1 - l/K)b{s) 

The Xi for the cut Pomerons are generated uniformly and then rescaled to 1 — x . Each 
cut Pomeron is then hadronized in its own center of mass using a modified independent 
fragmentation model with an energy dependent splitting function to reproduce the rise in 
dN/dy: 

f(x) = l-a + a(b(s) + l) b (s), b(s) = b + h log(s) 

The energy dependence is put into f(x) rather than P{K) because in the AKG scheme the 
single particle distribution comes only from the single chain. The probabilities for different 
flavors are taken to be 

u : d : s = .46 : .46 : .08 
to reproduce the experimental K/ir ratio. 
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3 Sample Jobs 

The simplest ISAJET job reads a user-supplied parameter file and writes a data file and a 
listing file. The following is an example of a parameter file which generates each type of 
event: 

SAMPLE TWOJET JOB 
800. ,100,2,50/ 
TWOJET 
PT 

50,100,50,100/ 
JETTYPE1 
'GL'/ 
JETTYPE2 

'UP' , 'UB' , 'DN' , 'DB' ,'ST> ,'SB'/ 
END 

SAMPLE DRELLYAN JOB 
800. ,100,2,50/ 
DRELLYAN 
QMW 

80,100/ 
WTYPE 
'W+' , 'W-'/ 
END 

SAMPLE MINBIAS JOB 
800. ,100,2,50/ 
MINBIAS 
END 

SAMPLE WPAIR JOB 
800. ,100,2,50/ 
WPAIR 
PT 

50,100,50,100/ 

JETTYPE1 

'W+' , 'W-' , 'Z0V 

JETTYPE2 

'W+' , 'W-' , 'Z0V 

WM0DE1 

'E+' , 'E-' , 'NUSV 

WM0DE2 

'QUARKS'/ 

END 

SAMPLE HIGGS JOB FOR SSC 

40000,100,1,1/ 

HIGGS 
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QMH 

400,1600/ 

HMASS 

800/ 

JETTYPE1 

'Z0'/ 

JETTYPE2 

'Z07 

WM0DE1 

'MU+' , 'MU-V 

WM0DE2 

'E+' , 'E-V 

PT 

50,20000,50,20000/ 
END 

SAMPLE SUSY JOB 

1800,100,1,10/ 

SUPERSYM 

PT 

50,100,50,100/ 
JETTYPE1 

'GLSS' , 'SQUARKSV 
JETTYPE2 

'GLSS' , ' SQUARKS ' / 
GAUGINO 
60,1,40,40/ 
S QUARK 

80.3,80.3,80.5,81.6,85,110/ 
FORCE 

29,30,1,-1/ 

FORCE 

21,29,1/ 

FORCE 

22,29,2/ 

FORCE 

23,29,3/ 

FORCE 

24,29,4/ 

FORCE 

25,29,5/ 

FORCE 

26,29,6/ 

END 

SAMPLE MSSM JOB FOR TEVATRON 



1800. ,100,1,1/ 

SUPERSYM 

BEAMS 

>P' , 'APV 

MSSMA 

200,-200,500,2/ 
MSSMB 

200,200,200,200,200/ 
MSSMC 

200,200,200,200,200,0,0,0/ 

JETTYPE1 

'GLSS'/ 

JETTYPE2 

'SQUARKS'/ 

PT 

100,300,100,300/ 
END 

SAMPLE MSSM SUGRA JOB FOR LHC 

14000,100,1,10/ 

SUPERSYM 

PT 

50,500,50,500/ 
SUGRA 

247,302,-617.5,10,-1/ 

TMASS 

175/ 

END 

SAMPLE SUGRA HIGGS JOB USING DEFAULT QMH RANGE 

14000,100,20,50/ 

HIGGS 

SUGRA 

200, 200,0, 2, +1/ 

HTYPE 

'HAO'/ 

JETTYPE1 

> GAUGINOS ' , ' SLEPTONS ' / 
JETTYPE2 

> GAUGINOS ' , ' SLEPTONS ' / 
END 

SAMPLE E+E- TO SUGRA JOB WITH POLARIZED BEAMS AND BREM/BEAMSTRAHLUNG 

500. ,100,1,1/ 

E+E- 

SUGRA 

125,125,0,3,1/ 
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TMASS 

175,-1,1/ 
EPOL 
-.9,0./ 
EEBEAM 

200. ,500. , .1072, .12/ 

JETTYPE1 

'ALL'/ 

JETTYPE2 

'ALL'/ 

NTRIES 

10000/ 

END 

SAMPLE WH JOB 

2000,100,0,0/ 

WHIGGS 

BEAMS 

'P' , 'AP'/ 

HMASS 

100./ 

JETTYPE1 

'W+' , 'W-\ 'HIGGS'/ 
JETTYPE2 

'W+' , 'W-' , 'HIGGS'/ 

WM0DE1 

'ALL'/ 

WM0DE2 

'ALL'/ 

PT 

10,300,10,300/ 
END 

SAMPLE EXTRA DIMENSIONS JOB 

14000,100,1,100/ 

EXTRADIM 

QMW 

5,1000/ 
QTW 

500,1000/ 
EXTRAD 

2,1000, .FALSE./ 
END 

SAMPLE ZJJ JOB AT LHC 

14000,100,1,100/ 

ZJJ 
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PT 

20 , 7000 , 20 , 7000 , 20 , 7000/ 

MIJLIM 

0,0,20,7000/ 

MTOT 

100,500/ 

NSIGMA 

200/ 

NTRIES 

10000/ 

END 

STOP 

See Section 6 of this manual for a complete list of the possible commands in a parameter 
file. Note that all input to ISAJET must be in UPPER case only. 

Subroutine RDTAPE is supplied to read events from an ISAJET data file, which is a 
machine-dependent binary file. It restores the event data to the FORTRAN common blocks 
described in Section 7. The skeleton of an analysis job using HBOOK and PAW from the 
CERN Program Library is provided in patch ISAPLT but is not otherwise supported. A 
Zebra output format based on code from the DO Collaboration is also provided in patch 
ISAZEB; see the separate documentation in patch ISZTEXT. 

3.1 DEC VMS 

On a VAX or ALPHA running VMS, ISAJET can be compiled by executing the .COM file 
contained in P=ISAUTIL,D=MAKEVAX. Extract this deck as ISAMAKE.COM and type 

@ISAMAKE 

This will run YPATCHY with the pilot patches described in Section 4 and the VAX flag to 
extract the source code, decay table, and documentation. The source code is compiled and 
made into a library, which is linked with the following main program, 

PROGRAM ISARUN 
C MAIN PROGRAM FOR ISAJET 

OPEN (UNIT=1,STATUS=' OLD' , F0RM= ' FORMATTED ' , READONLY) 

OPEN (UNIT=2 , STATUS= ' NEW ' , F0RM= ' UNFORMATTED ' ) 

OPEN (UNIT=3 , STATUS= ' OLD ' , F0RM= ' FORMATTED ' ) 

OPEN (UNIT=4 , STATUS= ' NEW ' , F0RM= ' FORMATTED ' ) 

CALL ISAJET (-1,2, 3, 4) 

STOP 

END 

to produce ISAJET.EXE. Two other executables, ISASUSY.EXE and ISASUGRA.EXE, 
will also be produced to calculate SUSY masses and decay modes without generating events. 
Temporary files can be removed by typing 
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OISAMAKE CLEAN 



Create an input file JOBNAME.PAR following the examples above or the instructions in 
Section 6 and run ISAJET with the command 

OISAJET JOBNAME 

using the ISAJET.COM file contained P=1SAUTIL,D=RUNVAX. This will create a binary 
output file JOBNAME . DAT and a listing file JOBNAME. LIS. Analyze the output data using the 
commands described in Section 8. 

There is also an simple interactive interface to ISAJET which will prompt the user for 
commands, write a parameter file, and optionally execute it. 

3.2 IBM VM/CMS 

On an IBM mainframe running VM/CMS, run YPATCHY with the pilot patches described 
in Section 4 and the IBM flag to extract the source code, decay table, and documentation. 
Compile the source code and link it with the main program 

PROGRAM ISARUN 
C MAIN PROGRAM FOR ISAJET 

OPEN (UNIT=1 , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 

OPEN (UNIT=2 , STATUS= ' NEW ' , FORM= ' UNFORMATTED ' ) 

OPEN (UNIT=3 , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 

OPEN (UNIT=4 , STATUS= ' NEW ' , FORM= ' FORMATTED ' ) 

CALL ISAJET (-1,2, 3, 4) 

STOP 

END 

to make ISAJET MODULE. 

Create a file called JOBNAME INPUT containing ISAJET input commands following the ex- 
amples above or the instructions in Section 6. Then run ISAJET using ISAJET EXEC, which 
is contained in P=ISAUTIL,D=RUNIBM. The events will be produced on JOBNAME DATA A 
and the listing on JOBNAME OUTPUT A. 

3.3 Unix 

The Makefile contained in P=ISAUTIL,D=MAKEUNIX has been tested on DEC Ultrix, 
Hewlett Packard HP-UX, IBM RS/6000 AIX, Linux, Silicon Graphics IRIX, Sun SunOS, 
and Sun Solaris. It should work with minor modifications on almost any Unix system with 
/bin/csh, ypatchy or nypatchy, and a reasonable Fortran 77 compiler. Extract the Makefile 
and edit it, changing the installation parameters to reflect your system. Note in particular 
that CERNlib is usually compiled with underscores postpended to all external names; you 
must choose the appropriate compiler option if you intend to link with it. Then type 

make 
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This should produce an executable isajet.x for the event generator, which links the code 
with the following main program: 

PROGRAM RUNJET 

CHARACTER*60 FNAME 

READ 1000, FNAME 
1000 FORMAT (A) 

PRINT 1020, FNAME 
1020 FORMAT (IX, 'Data file = ' ,A) 

OPEN (2 , FILE=FNAME , STATUS= ' NEW ' , FORM 

READ 1000, FNAME 

PRINT 1030, FNAME 
1030 FORMAT (IX, 'Parameter file = ' ,A) 

OPEN (3 , FILE=FNAME , STATUS= ' OLD ' , FORM 

READ 1000, FNAME 

PRINT 1040, FNAME 
1040 FORMAT (IX, 'Listing file = ' ,A) 

OPEN (4 , FILE=FNAME , STATUS= 'NEW , FORM 

READ 1000, FNAME 

OPEN ( 1 , FILE=FNAME , STATUS= ' OLD ' , FORM 
CALL ISAJET(-1,2,3,4) 
STOP 
END 

Two other executables, isasusy .x and isasugra.x, will also be produced to calculate SUSY 
masses and decay modes without generating events. Type 

make clean 

to delete the temporary files. 

Most Unix systems do not allow two jobs to read the same decay table file at the same 
time. The shell script in P=ISAUTIL,D=RUNUNIX copies the decay table to a temporary 
file to avoid this problem. Extract this file as isajet. Create an input file jobname.par 
following the examples above or the instructions in Section 6 and run ISAJET with the 
command 

isajet jobname 

This will create a binary output file jobname.dat and a listing file jobname. lis. Analyze 
the output data using the commands described in Section 8. 

This section only describes running ISAJET as a standalone program and generating 
output in machine-dependent binary form. The user may elect to analyze events as they are 
generated; this is discussed in Section 5 of this manual. 
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4 Patchy and PAM Organization 



Patchy is a code management system developed at CERN and used to maintain the CERN 
Library. It is used to provide versions of ISAJET for a wide variety of computers. Instructions 
for using PATCHY are available from http://wwwinfo.cern.ch/asdoc/Welcome.html. 

A master source file in Patchy is called a "PAM." The ISAJET PAM contains all the 
source code and documentation plus Patchy commands to include common blocks and to 
select the desired version. It is divided into the following patches: 

ISACDE: contains all common blocks, etc. These are divided into decks based on their 
usage. 

ISADATA: contains block data ALDATA. This must always be loaded when using ISAJET. 
ISAJET: contains the code for generating events. Each subroutine is in a separate deck 
with the same name. 

ISASSRUN: contains the main program for IS ASUS Y, which prompts for input parameters 
and prints out all the decay modes. It is selected by * ISASUSY. 

ISASUSY: contains code to calculate all the decay widths and branching fractions in the 
minimal supersymmetric model. 

I SAT APE: contains the code for reading and writing tapes, again with each subroutine on 
a separate deck. 

ISARUN: contains a main program and a simple interactive interface. It is selected by 
IF=INTERACT. 

ISAZEB: contains Zebra format output routines, an alternative to the ISATAPE routines. 
ISZRUN: contains the analog of ISAPLT for the Zebra format. 

ISAPLT: contains a simple calorimeter simulation and the skeleton of a histogramming 
job using HBOOK. 

ISATEXT: contains the instructions for using ISAJET, i.e. the text of this document. It 
also includes the documentation for ISASUSY. 

ISZTEXT: contains the instructions for the Zebra output routines and a description of the 
Zebra banks. 

ISADECAY: contains the input decay table. 

The code is actually maintained using RCS on a Silicon Graphics computer at BNL. 
Patchy is used primarily to handle common blocks and machine dependent code. 

The input to YPATCHY must contain both +USE cards, which define the wanted program 
version, and +EXE cards, which determine which patches or decks are written to the ASM 
file. To facilitate this selection, the ISAJET PAM contains the following pilot patches: 

* ISADECAY: USE selects ISADECAY and all corrections to it. 

*ISAJET: USE selects ISACDE, ISADATA, ISAJET, ISATAPE, ISARUN and all correc- 
tions to them. Note that ISARUN is not actually selected without +USE, INTERACT. 

*ISAPLT: USE selects ISACDE, ISAPLT, and all corrections to them. 

*ISASUSY: USE selects CDESUSY, ISASUSY, and ISASSRUN to create a program to 
calculate all the MSSM decay modes. 

*ISATEXT: USE selects ISACDE, ISATEXT, and all corrections to them. 

*ISAZEB: USE selects ISAJET with a Zebra output format. 

*ISZRUN: USE selects the Zebra analysis package. 
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Patches are provided to select the machine dependent features for specific computers or 
operating systems: 

ANSI: ANSI standard Fortran (no time or date functions) 

APOLLO: APOLLO - only tested by CERN 

CDC: CDC 7600 and 60-bit CYBER (obsolete) 

CRAY: CRAY with UNICOS 

DECS: DEC Station with Ultrix 

ETA: ETA 10 running Unix System V (obsolete) 

HPUX: HP/9000 7xx running Unix System V 

IBM: IBM 370 and 30xx running VM/CMS 

IBMRT: IBM RS/6000 running AIX 3.x or 4.x 

LINUX: PC running Linux with f2c/gcc or g77 compiler 

SGI: Silicon Graphics running IRIX 

SUN: Sun Sparcstation running SUNOS or Solaris 

VAX: DEC VAX or Alpha running VMS 
These patches in turn select a variety of patches and IF flags, allowing one to select more 
specific features, as indicated below. (Replace & by + everywhere.) 



WATCH, ANSI. 

&USE, DOUBLE. 

&USE,STDI0. 

&USE,M0VEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE , RANFCALL . 

&USE , NORANLUX , IF=-RANLUX . 

&USE , NOCERN , IF=-CERN . 

&E0D 



GENERIC ANSI FORTRAN. 
DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO RANLUX RANDOM NUMBERS. 
NO CERN LIBRARY. 



&PATCH, APOLLO . 

&DECK , BLANKDEK . 

&USE, DOUBLE. 

&USE,STDI0. 

&USE,M0VEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NORANLUX , IF=-RANLUX . 

&USE , NOCERN , IF=-CERN . 

&USE,IMPN0NE. 

&E0D. 



DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO RANLUX RANDOM NUMBERS. 
NO CERN LIBRARY. 
IMPLICIT NONE 



&PATCH,CDC. 
&USE, SINGLE. 
&USE,LEVEL2. 
&USE,CDCPACK. 
&USE, RANFCALL. 



CDC 7600 OR CYBER 175. 
SINGLE PRECISION. 
LEVEL 2 STORAGE. 

PACK 2 WORDS PER WORD FOR INPUT/OUTPUT. 
STANDARD RANSET AND RANGET CALLS. 
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&USE , NOCERN , IF=-CERN . 
&EOD 

&PATCH , CRAY . 

&USE, SINGLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , NOCERN , IF=-CERN . 

&EOD 

&PATCH,DECS. 

&USE,SUN. 

&EOD 

&PATCH , ETA . 

&USE, SINGLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFCALL . 

&USE , NOCERN , IF=-CERN . 

&EOD 

&PATCH,HPUX. 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NORANLUX , IF=-RANLUX . 

&USE , NOCERN , IF=-CERN . 

&USE,IMPNONE. 

&EOD 

&PATCH,IBM. 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NORANLUX , IF=-RANLUX . 

&USE , NOCERN , IF=-CERN . 

&EOD 

&PATCH, IBMRT . 



NO CERN LIBRARY. 



CRAY XMP OR 2. 
SINGLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
NO CERN LIBRARY. 



DEC STATION (ULTRIX) 



ETA- 10. 

SINGLE PRECISION. 
STANDARD FORTRAN 77 
FORTRAN REPLACEMENT 
STANDARD RANSET AND 
NO CERN LIBRARY. 



TAPE INPUT/ OUTPUT. 
FOR MOVLEV. 
RANGET CALLS. 



HP/9000 7XX RUNNING UNIX. 
DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO RANLUX RANDOM NUMBERS. 
NO CERN LIBRARY. 
IMPLICIT NONE 



IBM 370 OR 30XX. 
DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO RANLUX RANDOM NUMBERS. 
NO CERN LIBRARY. 



IBM RS/6000 WITH AIX 3.1 
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&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE , RANFCALL . 

&USE , NORANLUX , IF=-RANLUX . 

&USE , NOCERN , IF=-CERN . 

&USE,IMPNONE. 

&EOD 

&PATCH, LINUX. 
&USE, DOUBLE. 

&USE,STDIO. STANDARD FORTRAN 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NORANLUX , IF=-RANLUX . 

&USE , NOCERN , IF=-CERN . 

&USE,IMPNONE. 

&EOD 

&PATCH,SGI. 

SILICON GRAPHICS 4D/XX. 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NORANLUX , IF=-RANLUX . 

&USE , NOCERN , IF=-CERN . 

&EOD 

&PATCH,SUN. 

&USE, DOUBLE. 

&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE, RANFCALL. 

&USE , NORANLUX , IF=-RANLUX . 

&USE , NOCERN , IF=-CERN . 

&EOD 

&PATCH,VAX. 
&USE, DOUBLE. 



DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO RANLUX RANDOM NUMBERS. 
NO CERN LIBRARY. 
IMPLICIT NONE 



IBM PC WITH LINUX l.X 
DOUBLE PRECISION. 
77 TAPE INPUT/OUTPUT. 

FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO RANLUX RANDOM NUMBERS. 
NO CERN LIBRARY. 
IMPLICIT NONE 



DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO RANLUX RANDOM NUMBERS. 
NO CERN LIBRARY. 



SUN (SPARC) 
DOUBLE PRECISION. 

STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO RANLUX RANDOM NUMBERS. 
NO CERN LIBRARY. 



DEC VAX 11/780 OR 8600. 
DOUBLE PRECISION. 
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&USE,STDIO. 

&USE,MOVEFTN. 

&USE , RANFFTN , IF=-CERN . 

&USE , RANFCALL . 

&USE , NORANLUX , IF=-RANLUX . 

&USE , NOCERN , IF=-CERN . 

&USE,IMPNONE. 

&EOD 



STANDARD FORTRAN 77 TAPE INPUT/OUTPUT. 
FORTRAN REPLACEMENT FOR MOVLEV. 
FORTRAN RANF. 

STANDARD RANSET AND RANGET CALLS. 
NO RANLUX RANDOM NUMBERS. 
NO CERN LIBRARY. 
IMPLICIT NONE 



An empty patch INTERACT selects a main program and an interactive interface which 
will prompt the user for parameters and do some error checking. A patch CERN allows 
ISAJET to take the random number generator RANF and several other routines from the 
CERN Library; to use this include the Patchy command 

&USE,CERN. 

Similarly, a patch PDFLIB enables the interface to the PDFLIB parton distribution compi- 
lation by H. Plothow-Besch: 

&USE, PDFLIB 

The only internal links with PDFLIB are calls to the routines PDFSET, PFTOPDG, and 
DXPDF, and the common blocks W50510 and W50517, 

C Copy of PDFLIB common block 

C0MM0N/W50510/IFLPRT 

INTEGER IFLPRT 

SAVE /W50510/ 
C Copy of PDFLIB common block 

C0MM0N/W50517/N6 

INTEGER N6 

SAVE /W50517/ 

which are used to specify the level of output messages and the logical unit number for them. 
A patch RANLUX selects the RANLUX random number generator instead of the 48-bit 
congruential generator RANF: 

&USE , RANLUX 

RANLUX has a very long period, and each 32-bit integer seed gives an independent sequence. 

In general it should be sufficient to run YPATCHY with the following cradle (replace & 
with + everywhere): 



&USE , ( * I SA JET , * I S ATEXT , * I S ADECA Y , * I S APLT ) 
&USE, ANSI, DECS, HPUX, IBM, IBMRT, SGI, SUN, . . . 
& [USE, INTERACT] . 
& [USE, CERN.] 
&[USE,HB00K3.] 



CHOOSE ONE. 

CHOOSE ONE. 

FOR INTERACTIVE MODE. 

FOR CERN LIBRARY. 

HBOOK 3 FOR ISAPLT. 
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&EXE. 
&PAM. 
&QUIT. 

The input to YPATCHY can also contain changes by the user. It is suggested that these 
not be made permanent parts of the PAM to avoid possible conflicts with later corrections. 
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5 Main Program 

A main program for ISAJET is not supplied in the isajet . car file. To generate events and 
write them to disk, the user should provide a main program which opens the files and then 
calls subroutine ISAJET. In the following sample, i,j,m,n are arbitrary unit numbers. 
Main program for VMS: 

PROGRAM RUNJET 

C 

C MAIN PROGRAM FOR ISAJET ON BNL VAX CLUSTER. 

C 

0PEN(UNIT=i,FILE='$2$DUA14: [ISAJET . ISALIBRARY] DECAY . DAT ' , 
$STATUS= ' OLD ' , F0RM= ' FORMATTED ' , READONLY) 
OPEN (UNIT= j , FILE= ' my j ob . dat ' , STATUS= ' NEW ' , F0RM= ' UNFORMATTED ' ) 
OPEN (UNIT=m , FILE= ' my j ob . par ' , STATUS= 'OLD' , F0RM= ' FORMATTED ' ) 
OPEN (UNIT=n , FILE= ' my j ob . 1 i s ' , STATUS= ' NEW ' , F0RM= ' FORMATTED ' ) 

C 

CALL ISAJET(+-i,+-j ,m,n) 

C 

STOP 
END 

Main program for IBM (VM/CMS) 
PROGRAM RUNJET 

C 

C MAIN PROGRAM FOR ISAJET ON IBM ASSUMING FILES HAVE BEEN 

C OPENED WITH FILEDEF. 

C 

CALL ISAJET (+-i ,+- j ,m,n) 

C 

STOP 
END 

Main program for Unix; this is created by the standard Makefile: 
PROGRAM RUNJET 

C 

C Main program for ISAJET on Unix 

C 

CHARACTER*60 FNAME 

C 

C Open user files 

READ 1000, FNAME 
1000 FORMAT (A) 

PRINT 1020, FNAME 
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1020 FORMAT (IX, 'Data file = ' ,A) 

OPEN (2 , FILE=FNAME , STATUS= ' NEW ' , FORM= ' UNFORMATTED ' ) 

READ 1000, FNAME 

PRINT 1030, FNAME 
1030 FORMAT (IX, 'Parameter file = ' ,A) 

OPEN (3 , FILE=FNAME , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 

READ 1000, FNAME 

PRINT 1040, FNAME 
1040 FORMAT (IX, 'Listing file = ' ,A) 

OPEN (4 , FILE=FNAME , STATUS= 'NEW, FORM= ' FORMATTED ' ) 
C Open decay table 

READ 1000, FNAME 

OPEN ( 1 , FILE=FNAME , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 

C 

C Run ISA JET 

CALL ISAJET(-1,2,3,4) 

C 

STOP 
END 

The arguments of ISAJET are tape numbers for files, all of which should be opened by 
the main program. 

TAPEi: Decay table (formatted). A positive sign prints the decay table on the output 
listing. A negative sign suppress printing of the decay table. 

TAPEj: Output file for events (unformatted). A positive sign writes out both resonances 
and stable particles. A negative sign writes out only stable particles. 

TAPEm: Commands as defined in Section 6 (formatted). 

TAPEn: Output listing (formatted). 
In the sample jobs in Section 3, TAPEm is the default Fortran input, and TAPEn is the 
default Fortran output. 

5.1 Interactive Interface 

To use the interactive interface, replace the call to ISAJET in the above main program by 

CALL ISASET (+-i, +- j ,m,n) 
CALL ISAJET(+-i,+-j ,m,n) 

ISASET calls DIALOG, which prompts the user for possible commands, does a limited 
amount of error checking, and writes a command file on TAPEm. This command file is 
rewound for execution by ISAJET. A main program is included in patch ISARUN to open 
the necessary files and to call ISASET and ISAJET. 
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5.2 User Control of Event Loop 



If the user wishes to integrate ISAJET with another program and have control over the event 
generation, he can call the driving subroutines himself. The driving subroutines are: 

ISAINI(+-i,+-j ,m,n): initialize ISAJET. The arguments are the same as for subroutine 
ISAJET. 

ISABEG(IFL): begin a run. IFL is a return flag: IFL=0 for a good set of commands; 
IFL=1001 for a STOP; any other value means an error. 

ISAEVT(I , OK, DONE) generate event I. Logical flag OK signifies a good event (almost 
always .TRUE.); logical flag DONE signifies the end of a run. 

ISAEND: end a run. 

There are also subroutines provided to write standard ISAJET records, or Zebra records if 
the Zebra option is selected: 

ISAWBG to write a begin-of-run record, should be called immediately after ISABEG 
ISAWEV to write an event record, should be called immediately after ISAEVT 
ISAWND to write an end-of-run record, should be called immediately after ISAEND 
The control of the event loop is somewhat complicated to accomodate multiple evolu- 
tion and fragmentation as described in Section 11. Note in particular that after calling 
ISAEVT one should process or write out the event only if OK=.TRUE. The check on the 
DONE flag is essential if one is doing multiple evolution and fragmentation. The following 
example indicates how events might be generated, analyzed, and discarded (replace & by + 
everywhere) : 

PROGRAM SAMPLE 

C 

&SELF,IF=IMPNONE 

IMPLICIT NONE 

&SELF 

&CDE,ITAPES 

&CDE,IDRUN 

&CDE,PRIMAR 

&CDE,ISL00P 

C 

INTEGER JTDKY , JTEVT , JTCOM , JTLIS , IFL , IL00P 
LOGICAL OK, DONE 
SAVE IL00P 

c 

C> Open files as above 

C> Call user initialization 

c 

C 

C Initialize ISAJET 

C 

CALL ISAINI(-i,0,m,n) 
1 IFL=0 
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CALL ISABEG(IFL) 
IF(IFL.NE.O) STOP 

C 

C Event loop 

C 

IL00P=0 
101 CONTINUE 

IL00P=IL00P+1 
C Generate one event - discard if .NOT. OK 

CALL ISAEVT(IL00P,0K,D0NE) 
IF (OK) THEN 

c 

C> Call user analysis for event 

c 

END IF 

IF (.NOT. DONE) GO TO 101 

C 

C Calculate cross section and luminosity 

C 

CALL ISAEND 

c 

C> Call user summary 

c 

GO TO 1 
END 

5.3 Multiple Event Streams 

It may be desirable to generate several different kinds of events simultaneously to study 
pileup effects. While normally one would want to superimpose minimum bias or low-pt jet 
events on a signal of interest, other combinations might also be interesting. It would be very 
inefficient to reinitialize IS A JET for each event. Therefore, a pair of subroutines is provided 
to save and restore the context, i.e. all of the initialization information, in an array. The 
syntax is 

CALL CTX0UT(NC,VC,MC) 
CALL CTXIN(NC,VC,MC) 

where VC is a real array of dimension MC and NC is the number of words used, about 20000 
in the standard case. If NC exceeds MC, a warning is printed and the job is terminated. The 
use of these routines is illustrated in the following example, which opens the files with names 
read from the standard input and then superimposes on each event of the signal sample 
three events of a pileup sample. It is assumed that a large number of events is specified in 
the parameter file for the pileup sample so that it does not terminate. 

PROGRAM SAMPLE 
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c 

C Example of generating two kinds of events. 

C 

CHARACTER*60 FNAME 
REAL VC1 (20000) ,VC2 (20000) 
LOGICAL 0K1,D0NE1,0K2,D0NE2 
INTEGER NCI , NC2 , IFL , IL00P ,12, IL00P2 

C 

C Open decay table 

READ 1000, FNAME 
1000 FORMAT (A) 

OPEN ( 1 , FILE=FNAME , STATUS= 'OLD' , F0RM= ' FORMATTED ' ) 
C Open user files 

READ 1000, FNAME 

OPEN (3 , FILE=FNAME , STATUS= 'OLD' , F0RM= ' FORMATTED ' ) 
READ 1000, FNAME 

OPEN (4 , FILE=FNAME , STATUS= 'NEW , FORM= ' FORMATTED ' ) 
READ 1000, FNAME 

OPEN ( 13 , FILE=FNAME , STATUS= ' OLD ' , FORM= ' FORMATTED ' ) 
READ 1000, FNAME 

OPEN ( 14 , FILE=FNAME , STATUS= 'NEW , FORM= ' FORMATTED ' ) 

C 

C Initialize ISAJET 

CALL ISAINK-1,0,3,4) 
CALL CTXOUT (NCI, VC1, 20000) 
CALL ISAINK-1,0,13,14) 
IFL=0 

CALL ISABEG(IFL) 

IF(IFL.NE.O) ST0P1 

CALL CTXOUT (NC2,VC2, 20000) 

IL00P2=0 

CALL user_initialization_routine 

C 

1 IFL=0 

CALL CTXIN(NC1,VC1, 20000) 

CALL ISABEG(IFL) 

CALL CTXOUT (NCl.VC 1,20000) 

IF(IFL.NE.O) GO TO 999 

ILOOP=0 

C 

C Main event 

C 

101 CONTINUE 

IL00P=IL00P+1 
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CALL CTXIN (NCI ,VC1, 20000) 
CALL ISAEVT(IL00P,0K1,D0NE1) 
CALL CTXOUT (NCI, VC1, 20000) 
IF(.NOT.OKl) GO TO 101 
CALL user_analysis_routine 

C 

C Pileup 
C 

CALL CTXIN (NC2,VC2, 20000) 
12=0 

201 CONTINUE 

IL00P2=IL00P2+1 

CALL ISAEVT(IL00P2,0K2,D0NE2) 

IFC0K2) 12=12+1 

IF(D0NE2) ST0P2 

CALL user_analysis_routine 
IFCI2.LT. 3) GO TO 201 
CALL CTXOUT (NC2,VC2, 20000) 

C 

IF(.NOT.DONEl) GO TO 101 

C 

C Calculate cross section and luminosity 

C 

CALL CTXIN (NCI ,VC1 ,20000) 
CALL ISAEND 
GO TO 1 

C 

999 CALL CTXIN (NC2,VC2, 20000) 
CALL ISAEND 

CALL user_termination_routine 

STOP 

END 

It is possible to superimpose arbitrary combinations of events, including events of the same 
reaction type with different parameters. In general the number of events would be selected 
randomly based on the cross sections and the luminosity. 

At this time CTXOUT and CTXIN cannot be used with the Zebra output routines. 

5.4 Main Programs for ISASUSY and ISASUGRA 

In addition to the event generator, isajet.car contains two programs to calculate SUSY 
masses and decay modes: ISASUSY, which accepts weak scale parameters, and ISASUGRA, 
which calculates the weak scale parameters from those at some high scale. The main pro- 
grams, SSRUN and SUGRUN respectively, are included. They both prompt for interactive 
input and then call the appropriate ISAJET subroutines. The output is formatted and 



38 



printed by SSPRT and SUGPRT respectively. Executables for ISASUSY and ISASUGRA 
are built by the Unix Makefile and VMS isamake.com. 

It is fairly straightforward to modify these routines to scan SUSY parameters, but given 
the variety of possible scans, no attempt has been made to provide code for this. If only 
masses are needed, SSMSSM can be modified to remove the calls to the routines that calculate 
branching ratios. 
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6 Input 



6.1 Input Format 

ISAJET is controlled by commands read from the specified input file by subroutine READ IN. 
(In the interactive version, this file is first created by subroutine DIALOG.) Syntax errors 
will generate a message and stop execution. Based on these commands, subroutine LOGIC 
will setup limits for all variables and check for inconsistencies. Several runs with different 
parameters can be combined into one job. The required input format is: 

Title 

Ecm , Nevent , Npr int , N j ump/ 
Reaction 

(Optional parameters) 
END 

(Optional additional runs) 
STOP 

with all lines starting in column 1 and typed in upper case. These lines are explained below. 

Title line: Up to 80 characters long. If the first four letters are STOP, control is returned 
to main program. If the first four letters are SAME, the parameters from previous run are 
used excepting those which are explicitly changed. 

Ecm line: This line must always be given even if the title is SAME. It must give the 
center of mass energy (Ecm) and the number of events (Nevent) to be generated. One may 
also specify the number of events to be printed (Nprint) and the increment (Njump) for 
printing. The first event is always printed if Nprint > 0. For example: 

800. ,1000,10,100/ 

generates 1000 events at E cm = 800 GeV and prints 10 events. The events printed are: 

1,100,200, Note that an event typically takes several pages of output. This line is read 

with a list directed format (READ*). 

After Nprint events have been printed, a single line containing the run number, the event 
number, and the random number seed is printed every Njump events (if Njump is nonzero). 
This seed can be used to start a new job with the given event if in the new run NSIGMA is 
set equal to zero: 

SEED 
value/ 
NSIGMA 
0/ 

In general the same events will only be generated on the same type of computer. 

Reaction line: This line must be given unless title is SAME, when it must be omitted. 
It selects the type of events to be generated. The present version can generate TWOJET, 
E+E-, DRELLYAN, MINBIAS, WPAIR, SUPERSYM, HIGGS, PHOTON, TCOLOR, or 
WHIGGS events. This line is read with an A8 format. 



40 



6.2 Optional Parameters 



Each optional parameter requires two lines. The first line is a keyword specifying the param- 
eter and the second line gives the values for the parameter. The parameters can be given in 
any order. Numerical values are read with a list directed format (READ*), jet and particle 
types are read with a character format and must be enclosed in quotes, and logical flags with 
an LI format. All momenta are in GeV and all angles are in radians. 
The parameters can be classified in several groups: 



Jet Limits: 


W/H Limits: 


Decays: 


Constants: 


Other: 


JETTYPE1 


HTYPE 


FORCE 


AMSB 


BEAMS 


JETTYPE2 


PHIW 


FORCE1 


CUTJET 


EPOL 


JETTYPE3 


QMH 


NOB 


CUTOFF 


EEBEAM 


MIJLIM 


QMW 


NODECAY 


EXTRAD 


EEBREM 


MTOT 


QTW 


NOETA 


FRAGMENT 


NPOMERON 


P 


THW 


NOEVOLVE 


GAUGINO 


NSIGMA 


PHI 


WTYPE 


NOFRGMNT 


GMSB 


NTRIES 


PT 


XW 


NOGRAV 


GMSB2 


PDFLIB 


TH 


YW 


NOPIO 


HMASS 


SEED 


X 




NOTAU 


HMASSES 


STRUC 


Y 






LAMBDA 


WFUDGE 


WMODE1 






MGVTNO 


WMMODE 


WMODE2 






MSSMA 
MSSMB 


WPMODE 
ZOMODE 



MSSMC 

MSSMD 

MSSME 

NUSUG1 

NUSUG2 

NUSUG3 

NUSUG4 

NUSUG5 

SIGQT 

SIN2W 

SLEPTON 

SQUARE 

SSBCSC 

SUGRA 

SUGRHN 

TCMASS 

TMASS 

WMASS 



It may be helpful to know that the TWOJET, WPAIR, PHOTON, SUPERSYM, and 
WHIGGS processes use the same controlling routines and so share many of the same vari- 
ables. In particular, PT limits should normally be set for these processes, and JETTYPE1 
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and JETTYPE2 are used to select the reactions. Similarly, the DRELLYAN, HIGGS, and 
TCOLOR processes use the same controlling routines since they all generate s-channel reso- 
nances. The mass limits for these processes are set by QMW. Normally the QMW limits will 
surround the W ± , Z°, or Higgs mass, but this is not required. (QMH acts like QMW for the 
Higgs process.) For historical reasons, JETTYPE1 and JETTYPE2 are used to select the W 
decay modes in DRELLYAN, while WMODE1 and WMODE2 select the W decay modes for 
WPAIR, HIGGS, and WHIGGS. Also, QTW can be used to generate DRELLYAN events 
with non-zero transverse momentum, whereas HIGGS automatically fixes QTW to be zero. 
(Of course, non-zero transverse momentum will be generated by gluon radiation.) 
For example the lines 

P 

40. ,50. ,10. ,100./ 

would set limits for the momentum of jet 1 between 40 and 50 GeV, and for jet 2 between 
10 and 100 GeV. As another example the lines 

WTYPE 
'W+V 

would specify that for DRELLYAN events only W+ events will be generated. If for a 
kinematic variable only the lower limit is specified then that parameter is fixed to the given 
value. Thus the lines 

P 

40. ,,10./ 

will fix the momentum for jet 1 to be 40 GeV and for jet 2 to be 10 GeV. If only the upper 
limit is specified then the default value is used for the lower limit. Jet 1 or jet 2 parameters 
for DRELLYAN events refer to the W decay products and cannot be fixed. If QTW is fixed 
to 0, then standard Drell-Yan events are generated. 

A complete list of keywords and their default values follows. 
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Keyword 
Values 



Explanation 



Default values 



AMSB 

m ,m 3 /2,tan/3,sgn// 



none 



Anomaly-mediated SUSY breaking 
scalar mass, gravitino mass, 
VEV ratio, sign 



BEAMS 
typei,type 2 



'PVP' 



Initial beams. Allowed are 
'P','AP','N','AN'. 



CUTJET 

He 

CUTOFF 

/i 2 , v 

EEBEAM 

V^mm; ^f^rnax) T, <7 Z none 



6. 



.200,1.0 



EEBREM 
EPOL 

P L (e-),P L (e+) 



none 



0,0 



Cutoff mass for QCD jet 
evolution. 

Cutoff qt 2 = n 2 Q v for 
DRELLYAN events. 

impose brem/beamstrahlung 
min and max subprocess energy, 
beamstrahlung parameter T 
longitudinal beam size a z in mm 

impose bremsstrahlung for e + e~ 
min and max subprocess energy 

Polarization of e~ (e + ) beam, 
P L (e) = (n L - n R )/(n L - n R ), 
so that -1 < P L < 1 



EXTRAD 

<5,M D ,UVCUT 



None 



Parameters for EXTRADIM process 
UVCUT is logical flag 



FORCE 

Mi, -»»5/ 



FORCE1 

M'i, —,hl 



None 



None 



Force decay of particles, 
±i ±(il + ... + i5). 
Can call 20 times. 
See note for i = quark. 

Force decay i — > il + ... + i5. 
Can call 40 times. 
See note for i = quark. 



FRAGMENT 

Pud)- ■ ■ 

GAUGINO 

mi,m 2 ,m 3 , m 4 



.4, 



Fragmentation parameters. 
See also SIGQT, etc. 



Masses for g, 7, 
50,0,100,100 W + , and Z° 
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GMSB 

A m ,M m ,iV 5 
tan/3,sgn/x,C gr 



none 



GMSB2 

ft,5M 2 H JM 2 Hu ,D Y (M) 1,0,0,0 
N 5l ,N 52 ,N 53 N 5 



HMASS 

m 



GMSB messenger SUSY breaking, 
mass, number of 5 + 5, VEV 
ratio, sign, gravitino scale 

non-minimal GMSB parameters 
gaugino mass multiplier 
Higgs mass shifts, D-term mass 2 
indep. gauge group messengers 

Mass for standard Higgs. 



HMASSES 
mi,. . . ,m 9 

HTYPE 
'HLO'/ or... 

JETTYPE1 

'GL','UP',... 

JETTYPE2 
'GLVUP',... 

JETTYPE3 
'GLVUP',... 



LAMBDA 
A 



0,...,0 



none 



'ALL' 



'ALL' 



'ALL' 



Higgs meson masses for 
charges 0,0,0,0,0,1,1,2,2. 

One MSSM Higgs type ('HLO', 
'HH0', or 'HA0') 

) Select types for jets: 

)'ALL'; 'GL'; 'QUARKS'='UP', 

)'UB','DN','DB','ST','SB', 

)'CH','CB','BT','BB','TP', 

)'TB','X','XB','Y','YB'; 

)'LEPTONS'='E-','E+','MU-', 

)'MU+','TAU-','TAU+'; 'NUS'; 

)'GM','W+','W-','Z0' 

) See note for SUSY types. 

QCD scale 



MGVTNO 

-^gravitino 

MIJLIM 

^ i j i -^min i ^max 

MSSMA 
m(g),fi, 
m(A), tan (3 

MSSMB 

m(qi) ,m(d r ) ,m(u r 
m(h),m(e r ) 



10 20 GeV 



0,0,1 GeV,lGeV 



Required 



Required 



Gravitino mass - ignored for 
GMSB model 

Multimet mass limits 



MSSM parameters - 
Gluino mass, //, A mass, 
tan (3 

MSSM 1st generation - 

Left and right soft squark and 

slepton masses 
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MSSMC 

m(q 3 ) ,m(b r ) ,m(t r ) , Required 

m(/ 3 ),m(r r ), 

At,Ai>,A T 

MSSMD 

m(g 2 ) ,m>{sr) ,">Ti(cr) , from MSSMB 
m(l 2 ),m(m?i r ) 

MSSME 
Mi,M 2 

MTOT 

M min ,M max None 
NOB 

TRUE or FALSE FALSE 
NODECAY 

TRUE or FALSE FALSE 
NOETA 

TRUE or FALSE FALSE 
NOEVOLVE 

TRUE or FALSE FALSE 
NOGRAV 

TRUE or FALSE FALSE 
NOHADRON 

TRUE or FALSE FALSE 
NONUNU 

TRUE or FALSE FALSE 
NOPIO 

TRUE or FALSE FALSE 
NOTAU 

TRUE or FALSE FALSE 

NPOMERON 
n u n 2 1,20 



MSSM 3rd generation - 
Soft squark masses, slepton 
masses, and squark and slepton 
mixings 

MSSM 2nd generation - 

Left and right soft squark and 

slepton masses 



Mass range for multiparton 
processes 

Suppress B decays to use 
external package. 

Suppress all decays. 

Suppress eta decays. 

Suppress QCD evolution and 
hadronization. 

Suppress gravitino decays in 
GMSB model 

Suppress hadronization of 
jets and beam jets. 

Suppress Z° neutrino decays. 

Suppress 7r° decays. 

Suppress tau decays to use 
external package. 

Allow rii < n < 112 cut pomerons. 
Controls beam jet mult. 



MSSM gaugino masses - 
MSSMA + GUT Default is to scale from gluino 
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NSIGMA 

n 20 
NTRIES 

n 1000 
NUSUG1 

Mi,M 2 ,M 3 none 
NUSUG2 

A t ,Ab,A T none 
NUSUG3 

M Hd ,M Hu none 
NUSUG4 

M UL ,M dR ,M UR , none 
NUSUG5 

M h .M none 

M TL ,Mr R 

P 

(l),...,Pmax(3) l.,0.5#cm 

PDFLIB 

'namei',vali,. . . None 
PHI 

min (l),- • • ,0max(3) 0,27T 

PHIW 

0min,0 max 

PT or PPERP 

Pt,min(l)v • • )Pt,max(3) .05E cm ,.2E cm 

QMH 
QMW 

^min j Q'max 



Generate n unevolved events 
for SIGF calculation. 

Stop if after n tries 
cannot find a good event. 

Optional non-universal SUGRA 
gaugino masses 

Optional non-universal SUGRA 
A terms 

Optional non-universal SUGRA 
Higgs masses 

Optional non-universal SUGRA 
lst/2nd generation masses 

Optional non-universal SUGRA 
3rd generation masses 

Momentum limits for jets. 

CERN PDFLIB parton distribution 
parameters. See PDFLIB manual. 

Phi limits for jets. 
Phi limits for W. 



p t limits for jets. 

Default for TWOJET only. 



Mass limits for Higgs. 
mE cm ,.2E cm Equivalent to QMW. 

Mass limits for W. 

.05E cmr 2,E cm 
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QTW q t limits for W. Fix q t = 

<?t,min,<?t,max .l,.025£ cm for standard Drell-Yan. 

SEED Seed <281474976710656.D0 for 

real/integer RANF or <2 31 for RANLUX. 

SIGQT Internal k t parameter for 

a .35 jet fragmentation. 

SIN2W Weinberg angle. See WMASS. 

sin 2 (^) .232 

SLEPTON Masses for z> e , e, z> M , jl, z> T , f 

m lr ..,rriQ 100,. .. ,101.8 

SQUARE Masses for u, d, s, c, b, t 

mi,. ..,m 6 100.3,. ..,240. 

SSBCSC Alternate mass scale for RGE 

M Mqut boundary conditions. 

STRUC Structure functions. CTEQ5L, 

name 'CTEQ5L' CTEQ3L, CTEQ2L, EHLQ, OR DO 



SUGRA 

m ,mi/ 2 ,Ao, none 
tan /5,sgn jj, 
TH or THETA 
6» min (l),. . . ,6» max (3) 0,tt 

SUGRHN 

m VT ,M N ,A n ,m i>R 0,1^20,0,0 



THW 

^min^max 0,7T 

TCMASS 

m,r 1000,100 
TMASS 

mt,m y ,m x 180.,-1.,-1. 



Minimal supergravity parameters 
scalar M, gaugino M, trilinear 
breaking term, vev ratio, +-1 
Theta limits for jets. Do not 
also set Y. 

SUGRA see-saw //-effect 
nu-mass, int. scale, 
GUT scale nu SSB terms 

Theta limits for W. Do not 
also set YW. 

Technicolor mass and width. 



t, y, and x quark masses. 
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WFUDGE 

factor 



1.85 



Fudge factor for DRELLYAN 
evolution scale. 



WMASS 

M W ,M Z 80.2, 91.19 

WMMODE 

'UP',. . . ,'TAU+' 'ALL' 
WMODE1 

'UP','UB',... 'ALL' 
WMODE2 

'UP','UB',... 'ALL' 
WPMODE 

'UP',. . . ,'TAU+' 'ALL' 
WTYPE 

type^typea 'GM','Z0' 
X 

•^min(l)v • • )-^max(3) 1,1 



W and Z masses. See SIN2W. 



Decay modes for W in parton 
cascade. See JETTYPE. 

) 

) Decay modes for WPAIR. 
)Same code for quarks and 
)leptons as JETTYPE. 

) 

Decay modes for W + in parton 
cascade. See JETTYPE. 

Select W type: W+,W-,GM,Z0. 
Do not mix W+,W- and GM,Z0. 

Feynman x limits for jets. 



XGEN Jet fragmentation, Peterson 

a(l),...,a(8) .96,3,0,.8,.5,... with e = a(n)/m 2 , n = 4-8. 

XGENSS Fragmentation of GLSS, UPSS, 

a(l),. . . ,a(7) .5, .5,... etc. with e = a(n)/m * *2 

XW Feynman x limits for W. 

•£min;3'max ljl 

Y Y limits for each jet. 

2/min(l)v • • ,2/max(3) from PT Do not also set TH. 

YW Y limits for W. 

ymin,ymax from QTW,QMW Do not set both YW and THW. 

Z0MODE Decay modes for Z° in parton 



'UP',. . . ,'TAU+' 'ALL' cascade. See JETTYPE. 
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6.3 Kinematic and Parton-type Parameters 

While the TWO JET PT limits and the DRELLYAN QMW limits are formally optional 
parameters, they are set by default to be fractions of y/s. Thus, for example, the parameter 
file 

DEFAULT TWO JET JOB 

14000,100,1,100/ 

TWO JET 

END 

STOP 

will execute, but it will generate jets between 5% and 20% of y/s, which is probably not 
what is wanted. Similarly, the parameter file 

DEFAULT DRELLYAN JOB 

14000,100,1,100/ 

DRELLYAN 

END 

STOP 

will generate 7 + Z events with masses between 5% and 20% of y/s, not masses around the 
Z mass, and transverse momenta between 1 GeV and 2.5% of y/s. 

Normally the user should set PT limits for TWOJET, PHOTON, WPAIR, SUPERSYM, 
and WHIGGS events and QMW and QTW limits for DRELLYAN, HIGGS, and TCOLOR 
events. If these limits are not set, they will be selected as fractions of E cm . This can give 
nonsense. For TWOJET the p t range should usually be less than about a factor of two 
except for b and t jets at low p t to produce uniform statistics. For W + , W~ , or Z° events 
or for Higgs events the QMW (QMH) range should usually include the mass. But one can 
select different limits to study, e.g., virtual W production or the effect of a lighter or heavier 
Higgs on WW scattering. If only t decays are selected, then the lower QMW limit must be 
above the t threshold. For standard Drell-Yan events QTW should be fixed to zero, 

QTW 

0/ 

Transverse momenta will then be generated by initial state gluon radiation. A range of 
QTW can also be given. For SUPERSYM either the masses and decay modes should be 
specified, or the MSSM, SUGRA, GMSB, or AMSB parameters should be given. For fourth 
generation quarks it is necessary to specify the quark masses. 

Note that if the limits given cover too large a kinematic range, the program can become 
very inefficient, since it makes a fit to the cross section over the specified range. NTRIES 
has to be increased if narrow limits are set for X, XW or for jet 1 and jet 2 parameters 
in DRELLYAN events. For larger ranges several runs can be combined together using the 
integrated cross section per event SIGF/NEVENT as the weight. This cross section is 
calculated for each run by Monte Carlo integration over the specified kinematic limits and 
is printed at the end of the run. It is corrected for JETTYPEi, WTYPE, and WMODEi 
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selections; it cannot be corrected for branching ratios of forced decays or for WPMODE, 
WMMODE, or ZOMODE selections, since these can affect an arbitrary number of particles. 

To generate events over a large range, it is much more efficient to combine several runs. 
This is facilitated by using the special job title SAME as described above. Note that 
SAME cannot be used to combine standard DRELLYAN events (QTW fixed equal to 0) 
and DRELLYAN events with nonzero QTW. 

The cross sections for multiparton final states in general have infrared and collinear 
singularities. To obtain sensible results, it is in general essential to set limits both on the pt 
of each final parton using PT and on the mass of each pair of partons using M1JL1M. The 
default lower limits are all 1 GeV. Using these default limits without thought will likely give 
absurd results. 

For TWOJET, DRELLYAN, and most other processes, the JETTYPEi and WTYPEi 
keywords should be used to select the subprocesses to be included. For e + e~ — > W + W~, 
Z°Z°, use FORCE and FORCE1 instead of WMODEi to select the W decay modes. Note 
that these do not change the calculated cross section. (In the E+E- process, the W and Z 
decays are currently treated as particle decays, whereas in the WPAIR and HIGGS processes 
they are treated as 2 — > 4 parton processes.) 

For HIGGS with W + W~ or Z°Z° decays allowed it is generally necessary to set PT limits 
for the W's, e.g. 

PT 

50,20000,50,20000/ 

If this is not done, then the default lower limit of 1 GeV is used, and the t-channel exchanges 
will dominate, as they should in the effective W approximation. Depending on the other 
parameters, the program may fail to generate an event in NTRIES tries. 

6.4 SUSY Parameters 

SUPERSYM (SUSY) by default generates just gluinos and squarks in pairs. There are no 
default masses or decay modes. Masses can be set using GAUGINO, SQUARE, SLEPTON, 
and HMASSES. Decay modes can be specified with FORCE or by modifying the decay 
table. Left and right squarks are distinguished but assumed to be degenerate, except for 
stops. Since version 7.11, types must be selected with JETTYPEi using the supersymmetric 
names, e.g. 

JETTYPEI 

'GLSS' , 'UPSSL' , 'UPSSRV 

Use of the corresponding standard model names, e.g. 

JETTYPEI 
'GL' , 'UP'/ 

and generation of pure photinos, winos, and zinos are no longer supported. 

If MSSMA, MSSMB and MSSMC are given, then the specified parameters are used to 
calculate all the masses and decay modes with the ISASUSY package assuming the minimal 
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supersymmetric extension of the standard model (MSSM). There are no default values, so 
you must specify values for each MSSMi, i=A-C. MSSMD can optionally be used to set the 
second generation squark and slepton parameters; if it is omitted, then the first generation 
ones are used. MSSME can optionally be used to set the U(l) and SU(2) gaugino masses; 
if it is omitted, then the grand unification values are used. The parameters and the use of 
the MSSM is preserved if the title is SAME. FORCE can be used to override the calculated 
branching ratios. 

The MSSM option also generates charginos and neutralinos with cross sections based on 
the MSSM mixing angles in addition to squarks and sleptons. These can be selected with 
JETTYPEi; the complete list of supersymmetric options is: 

'GLSS' , 

'UPSSL' , 'UBSSL' , 'DNSSL' , 'DBSSL' , 'STSSL' , 'SBSSL' , 'CHSSL' , 'CBSSL' , 
'BTSS1' , 'BBSS1' , 'TPSS1' , 'TBSS1' , 

'UPSSR' , 'UBSSR' , 'DNSSR' , 'DBSSR' , 'STSSR' , 'SBSSR' , 'CHSSR' , 'CBSSR' , 
'BTSS2' , 'BBSS2 ' , 'TPSS2 ' , 'TBSS2' , 

'W1SS+' , ' W1SS- ' , 'W2SS+' , 'W2SS-' , 'Z1SS' , 'Z2SS' , 'Z3SS' , 'Z4SS' , 

'NUEL' , ' ANUEL ' , >EL-> , >EL+' , ' NUML ' , ' ANUML ' ,MUL-' , ' MUL+ ' , ' NUTL ' , 

' ANUTL ' , 'TAU1- ' , 'TAU1+' , ' ER- ' , 'ER+' , 'MUR-' , 'MUR+' , 'TAU2-' , 'TAU2+' , 

'ZO' , 'HLO' , 'HHO' , 'HAO' , J H+' , 'H-' , 

' SQUARKS ' , ' GAUGINOS ' , ' SLEPTONS ' , ' ALL ' . 

Note that mixing between L and R stop states results in 1 (light) and 2 (heavy) stop, sbottom 
and stau eigenstates, which depend on the input parameters of left- and right- scalar masses, 
plus A terms, fx and tan/3. The last four JETTYPE's generate respectively all allowed 
combinations of squarks and antisquarks, all combinations of charginos and neutralinos, all 
combinations of sleptons and sneutrinos, and all SUSY particles. 

For SUSY Higgs pair production or associated production in E+E-, select the appropriate 
JETTYPE's, e.g. 

JETTYPEI 
'ZOV 
JETTYPE2 
'HLO'/ 

As usual, this gives only half the cross section. For single production of neutral SUSY 
Higgs in pp and pp reactions, use the HIGGS process together with the MSSMi, SUGRA, 
GMSB, or AMSB keywords. You must specify one and only one Higgs type using 

HTYPE 

'HLO' or 'HHO' or 'HAO'/ ««< One only! 

If no QMH range is given, one is calculated using M ± 5T for the selected Higgs. Decays 
into quarks, leptons, gauge bosons, lighter Higgs bosons, and SUSY particles are generated 
using the on-shell branching ratios from ISASUSY. You can use JETTYPEi to select the 
allowed Higgs modes and WMODEi to select the allowed decays of W and Z bosons. Since 
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heavy SUSY Higgs bosons couple weakly to W pairs, WW fusion and WW scattering are 
not included. 

SUGRA can be used instead of MSSMi to generate MSSM decays with parameters de- 
termined from mo, mi/2, Aq, tan/3, and sgn/i = ±1 in the minimal supergravity framework. 
The NUSUGi keywords can optionally be used to specify additional parameters for non- 
universal SUGRA models, while SUGRHN is used to specify the parameterf of an optional 
right-handed neutrino. Similarly, the GMSB keyword is used to specify the A, M m , N 5 , 
tan/5, sgn/i = ±1, and C grav parameters of the minimal Gauge Mediated SUSY Breaking 
model. GMSB2 can optionally be used to specify additional parameters of non-minimal 
GMSB models. The AMSB keyword is used to specify mo, m 3 / 2 , tan/3, and sgn/i for the 
minimal Anomaly Mediated SUSY Breaking model. Note that m 3 / 2 is much larger than the 
weak scale, typically 50 TeV. 

WHIGGS is used to generate W plus neutral Higgs events. For the Standard Model 
the JETTYPE is HIGGS. If any of the SUSY models is specified, then the appropriate SUSY 
Higgs type should be used, most likely HLO. In either case WMODEi is used to specify the W 
decay modes. The Higgs is treated as a particle; its decay modes can be set using FORCE. 

6.5 Forced Decay Modes 

The FORCE keyword requires special care. Its list must contain the numerical particle 
IDENT codes, e.g. 

FORCE 

140,130,-120/ 

The charge-conjugate mode is also forced for its antiparticle. Thus the above example forces 
both D° — > K + n~ and D° — > K~n + . If only a specific decay is wanted one should use the 
FORCE1 command; e.g. 

F0RCE1 

140,130,-120/ 

only forces D° — > K + n~ . 

To force a heavy quark decay one must generally separately force each hadron containing 
it. If the decay is into three leptons or quarks, then the real or virtual W propagator is 
inserted automatically. Since Version 7.30, top and fourth generation quarks are treated as 
particles and decayed directly rather than first being made into hadrons. Thus for example 

F0RCE1 
6,-12,11,5/ 

forces all top quarks to decay into an positron, neutrino and a b-quark (which will be 
hadronized). For the physical top mass, the positron and neutrino will come from a real W. 
Note that forcing t — > W + b and W + — > e + u e does not give the same result; the first uses the 
correct V — A matrix element, while the second decays the W according to phase space. 

Forced modes included in the decay table or generated by ISASUSY will automatically 
be put into the correct order and will use the correct matrix element. Modes not listed in 
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the decay table are allowed, but caution is advised because a wrong decay mode can cause 
an infinite loop or other unexpected effects. 

FORCE (FORCE1) can be called at most 20 (40) times in any run plus all subsequent 
'SAME' runs. If it is called more than once for a given parent, all calls are listed, and the 
last call is used. Note that FORCE applies to particles only, but that for gamma, W+, W-, 
Z0 and supersymmetric particles the same IDENT codes are used both as jet types and as 
particles. 

6.6 Parton Distributions 

The default parton distributions are fit CTEQ5L from the CTEQ Collaboration using lowest 
order QCD. The CTEQ3L, CTEQ2L, and the old EHLQ and Duke-Owens distributions can 
be selected using the STRUC keyword. 

If PDFLIB support is enabled (see Section 4), then any of the distributions in the PDFLIB 
compilation by H. Plothow-Besch can be selected using the PDFLIB keyword and giving the 
proper parameters, which are identical to those described in the PDFLIB manual and are 
simply passed to the routine PDFSET. For example, to select fit 29 (CTEQ3L) by the CTEQ 
group, leaving all other parameters with their default values, use 

PDFLIB 
'CTEQ' ,29D0/ 

Note that the fit-number and the other parameters are of type DOUBLE PRECISION 
(REAL on 64-bit machines). There is no internal passing of parameters except for those 
which control the printing of messages. 

6.7 Multiparton Processes 

For multiparton final states one should in general set limits on the total mass MTOT of the 
final state, on the minimum PT of each light parton, and on the minimum mass MIMLIM of 
each pair of light partons. Limits for PT are set in the ususal way. Limits for the mass 
of partons i,j are set using 

MIJLIM 

i , j ,Mmin,Mmax 

If i — j — 0, the limit is applied to all jet pairs. For example the following parameter file 
generates ZJJ events at the LHC with a mimimum p T of 20 GeV and a minimum mass of 
20 GeV for all jet pairs: 

GENERATE ZJJ with PTMIN = 20 GEV AND MMIN = 20 GEV 

14000,100,1,100/ 

ZJJ 

PT 

20 , 7000 , 20 , 7000 , 20 , 7000/ 
MIJLIM 



53 



0,0,20,7000/ 

MTOT 

100,500/ 

NSIGMA 

200/ 

NTRIES 

10000/ 

END 

STOP 

The default lower limits for PT and MIJLIM are 1 GeV. While these limits are sufficient to 
make the cross sections finite, they will in general not give physically sensible results. Thus, 

the user must think carefully about what limits should be set. 
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7 Output 



The output tape or file contains three types of records. A beginning record is written by 
a call to ISAWBG before generating a set of events; an event record is written by a call to 
ISAWEV for each event; and an end record is written for each run by a call to ISAWND. 
These subroutines load the common blocks described below into a single 

COMMON/ZEVEL/ZEVEL ( 1024) 

and write it out when it is full. A subroutine RDTAPE, described in the next section, inverts 
this process so that the user can analyze the event. 

ZEVEL is written out to TAPEj by a call to BUFOUT. For the CDC version IF = 
PAIRPAK is selected; BUFOUT first packs two words from ZEVEL into one word in 

C0MM0N/ZV0UT/ZV0UT(512) 

using subroutine PAIRPAK and then does a buffer out of ZVOUT to TAPEj. Typically at 
least two records are written per event. For all other computers IF=STDIO is selected, and 
ZEVEL is written out with a standard FORTRAN unformatted write. 



7.1 Beginning Record 

At the start of each run ISAWBG is called. It writes out the following common blocks: 

COMMON/DYLIM/QMIN , QMAX , QTMIN, QTMAX , YWMIN , YWMAX , XWMIN , XWMAX , THWMIN , 

2 THWMAX,PHWMIN,PHWMAX 

3 ,SETLMQ(12) 
SAVE /DYLIM/ 
LOGICAL SETLMQ 
EQUIVALENCE (BLIMl(l) ,QMIN) 

REAL QMIN , QMAX , QTMIN , QTMAX , YWMIN , YWMAX , XWMIN , XWMAX , THWMIN , 

+ THWMAX,PHWMIN,PHWMAX,BLIM1(12) 

QMIN, QMAX = W mass limits 

QTMIN,QTMAX = W q t limits 

YWMIN, YWMAX = W rj rapidity limits 

XWMIN,XWMAX = W x F limits 

THWMIN, THWMAX = W 6 limits 

PHWMIN,PHWMAX = W <p limits 

COMMON/ IDRUN/ IDVER, IDG (2) , IEVT, IEVGEN 
SAVE /IDRUN/ 

INTEGER IDVER , IDG , IEVT , IEVGEN 
IDVER = program version 

IDG(l) = run date (10000 xmonth+ 100 xday+year) 
IDG(2) = run time (10000 xhour+lOOxminute+second) 
IEVT = event number 
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1 Jet limits 

INTEGER MXLIM 
PARAMETER (MXLIM=8) 
INTEGER MXLX12 
PARAMETER (MXLX12=12*MXLIM) 

COMMON/ JETLIM/PMIN (MXLIM) ,PMAX (MXLIM) ,PTMIN (MXLIM) , PTMAX (MXLIM) , 
$YJMIN (MXLIM) ,YJMAX (MXLIM) , PHIMIN (MXLIM) ,PHIMAX (MXLIM) , 
$XJMIN (MXLIM) ,XJMAX (MXLIM) , THMIN (MXLIM) , THMAX (MXLIM) , 
$SETLMJ(12*MXLIM) 

SAVE /JETLIM/ 

COMMON/FIXPAR/FIXP (MXLIM) ,FIXPT (MXLIM) ,FIXYJ(MXLIM) , 
$FIXPHI (MXLIM) ,FIXXJ(MXLIM) ,FIXQM,FIXQT,FIXYW,FIXXW,FIXPHW 
SAVE /FIXPAR/ 

COMMON/SGNPAR/CTHS (2 , MXLIM) , THS (2 , MXLIM) , Y JS (2 , MXLIM) , X JS (2 , MXLIM) 
SAVE /SGNPAR/ 

REAL PMIN , PMAX , PTMIN , PTMAX , YJMIN , YJMAX , PHIMIN , PHIMAX , X JMIN , 

+ X JMAX , THMIN , THMAX , BLIMS ( 12*MXLIM) , CTHS , THS , Y JS , X JS 

LOGICAL SETLMJ 

LOGICAL FIXQM, FIXQT , FIXYW , FIXXW , FIXPHW 
LOGICAL FIXP , FIXPT , FIXY J , FIXPHI , FIXXJ 
EQUIVALENCE (BLIMS ( 1 ) , PMIN ( 1 ) ) 

PMIN,PMAX = jet momentum limits 

PTMIN, PTMAX = jet p t limits 

YJMIN,YJMAX = jet r) rapidity limits 

PHIMIN,PHIMAX = jet limits 

THMIN, THMAX = jet 9 limits 

INTEGER MXKEYS 

PARAMETER (MXKEYS=20) 

COMMON/KEYS/ IKEYS , KEYON , KEYS (MXKEYS) 

COMMON/XKEYS/REAC 

SAVE /KEYS/ , /XKEYS/ 

LOGICAL KEYS 

LOGICAL KEYON 

CHARACTER*8 REAC 

INTEGER IKEYS 
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KEYON = normally TRUE, FALSE if no good reaction 
KEYS = TRUE if reaction I is chosen 

1 for TWOJET 

2 for E+E- 

3 for DRELLYAN 

4 for MINBIAS 

5 for SUPERSYM 

6 for WPAIR 

REAC = character reaction code 

COMMON/PRIMAR/N JET , SCM , HALFE , ECM , IDIN (2) , NEVENT , NTRIES , NSIGMA 
SAVE /PRIMAR/ 

INTEGER N JET , IDIN , NEVENT , NTRIES , NSIGMA 
REAL SCM, HALFE, ECM 

NJET = number of jets per event 

SCM = square of com energy 

HALFE = beam energy 

ECM = com energy 

IDIN = ident code for initial beams 

NEVENT = number of events to be generated 

NTRIES = maximum number of tries for good jet parameters 

NSIGMA = number of extra events to determine SIGF 



INTEGER MXG0Q,MXG0J 
PARAMETER (MXG0Q=85,MXG0J=8) 

C0MM0N/Q1Q2/G0Q(MXG0Q,MXG0J) , GOALL (MXGO J) ,G0DY(4) , STDDY , 
$G0WW(25,2) ,ALLWW(2) , GOWMOD (25 ,MXG0 J) 
SAVE /Q1Q2/ 

LOGICAL GOQ , GOALL , GODY , STDDY , GOWW , ALLWW , GOWMOD 



GOQ(I,K) 



GOALL(K) 

GODY(I) 

1= 1 2 3 4 

GM W+ W- 

STDDY 

GOWW(I,K) 

ALLWW(K) 



ZO 



TRUE if quark type I allowed for jet k 
1=1234 5_6 789 10 11 12 13 

=^-guuddssccbbti 
I = 14 15 16 17 18 19 20 21 22 23 24 25 

=>• v e v e e" e + fT fi + v T v T r~ r" 

TRUE if all jet types allowed 
TRUE if W type I is allowed. 



TRUE if standard DRELLYAN 

TRUE if I is allowed in the decay of K for WPAIR. 

TRUE if all allowed in the decay of K for WPAIR. 



COMMON/QCDPAR/ALAM , ALAM2 , CUT JET , ISTRUC 
SAVE /QCDPAR/ 
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INTEGER ISTRUC 

REAL ALAM , ALAM2 , CUT JET 



QCD scale A 
QCD scale A 2 

cutoff for generating secondary partons 

3 for Eichten (EHLQ), 

4 for Duke (DO) 

5 for CTEQ 2L 

6 for CTEQ 3L 
-999 for PDFLIB 

COMMON/QLMASS/ AMLEP ( 100) , NQLEP , NMES , NBARY 
SAVE /QLMASS/ 

INTEGER NQLEP , NMES , NBARY 
REAL AMLEP 

AMLEP(6:8) = t,y,x masses, only elements written 



ALAM 
ALAM2 
CUTJET = 
ISTRUC = 



7.2 Event Record 

For each event ISAWEV is called. It writes out the following common blocks: 

COMMON/FINAL/NKINF , SIGF , ALUM , ACCEPT , NRECS 

SAVE /FINAL/ 

INTEGER NKINF, NRECS 

REAL SIGF, ALUM, ACCEPT 

SIGF = integrated cross section, only element written 

COMMON/ IDRUN/IDVER, IDG (2) , IEVT , IEVGEN 
SAVE /IDRUN/ 

INTEGER IDVER , IDG , IEVT , IEVGEN 

IDVER = program version 
IDG = run identification 
IEVT = event number 

COMMON/ JETPAR/P (3) ,PT(3) ,YJ(3) ,PHI(3) ,XJ(3) ,TH(3) , CTH (3) , STH(3) 

1 , JETTYP(3) ,SHAT,THAT,UHAT,QSQ,X1,X2,PBEAM(2) 

2 , QMW , QW , QTW , YW , XW , THW , QTMW , PHI W , SHAT 1 , THAT 1 , UHAT 1 , JWTYP 

3 ,ALFQSQ,CTHW,STHW,QOW 

4 ,INITYP(2) ,ISIGS,PBEAMS(5) 
SAVE /JETPAR/ 

INTEGER JETTYP , JWTYP , INITYP , ISIGS 

REAL P, PT,YJ, PHI, XJ,TH, CTH, STH, SHAT, THAT, UHAT, QSQ, XI, X2, 

+ PBEAM , QMW , QW , QTW , YW , XW , THW , QTMW , PHIW , SHAT1 , THAT1 , UHAT1 , 

+ ALFQSQ , CTHW , STHW , QOW , PBEAMS 
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1 ) 

r 


= jet momentum \p\ 


r 1 


= jet p t 


Y J 


= jet i] rapidity 


PHI 


= jet 


XJ 


= jet xf 


TH 


= jet 9 


CTH 


= jet cos(6>) 


STH 


= jet sin(#) 


JETTYP 


= jet type. The code is listed under /Q1Q2/ above 
continued. . . 



CITTAT TTT AT TTTTAT 




O T» r\ ri / * \~ \~ /~\T*1 n fT C* ~f~ HI 

IldlQ bCdXXt^IlIlg 6, 6, (/ 






nff'nr'f i ^ ro /O ^ 


X1 X9 




1111 Lildl JJdl liUll JU ]? 


PBEAM 




remaining beam momentum 


QMW 




W mass 


QW 




W momentum 


QTW 




W transverse momentum 


YW 




W rapidity 


XW 




W x F 


THW 




W 9 


QTMW 




Jfiw + Q 2 


PHIW 




W <f) 


SHAT1,THAT1,UHAT1 




invariants for W decay 


JWTYP 




W type. The code is listed under /Q1Q2/ above. 


ALFQSQ 




QCD coupling a s (Q 2 ) 


CTHW 




W cos(#) 


STHW 




W sm(9) 


QOW 




W energy 



INTEGER MXJSET,JPACK 

PARAMETER (MXJSET=400 , JPACK=1000) 

C0MM0N/JETSET/NJSET,PJSET(5,MXJSET) , JORIG(MXJSET) , JTYPE(MXJSET) , 
$JDCAY(MXJSET) 
SAVE /JETSET/ 

INTEGER NJSET , JORIG , JTYPE , JDCAY 
REAL PJSET 
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NJSET 

PJSET(1,I) 

PJSET(2,I) 

PJSET(3,I) 

PJSET(4,I) 

PJSET(5,I) 

JORIG(I) 



JTYPE(I) 
JDCAY(I) 

MXJSET 
JPACK 



number of partons 
p x of parton I 
p y of parton I 
p z of parton I 
Po of parton I 
mass of parton I 

JPACK* JET+K if I is a decay product of K. 
IF K=0 then I is a primary parton. 
(JET = 1,2,3 for final jets.) 
(JET = 11,12 for initial jets.) 
IDENT code for parton I 

JPACK*K1+K2 if Kl and K2 are decay products of I. 
If JDCAY(I)=0 then I is a final parton 
dimension for /JETSET/ arrays, 
packing integer for /JETSET/ arrays. 



INTEGER MXSIGS , IOPAK 

PARAMETER (MXSIGS=3000 , I0PAK=100) 

COMMON/ JETSIG/SIGMA , SIGS (MXSIGS) ,NSIGS , INOUT(MXSIGS) ,SIGEVT 

SAVE /JETSIG/ 

INTEGER NSIGS,INOUT 

REAL SIGMA, SIGS, SIGEVT 

SIGMA = cross section summed over types 

SIGS (I) = cross section for reaction I (not written) 

NSIGS = number of nonzero cross sections (not written) 

INOUT(I) = packed partons for process I (not written) 

MXSIGS = dimension for JETSIG arrays (not written) 

SIGEVT = partial cross section for selected channel 



INTEGER MXPTCL,IPACK 

PARAMETER (MXPTCL=4000 , IPACK=10000) 

C0MM0N/PARTCL/NPTCL,PPTCL(5,MXPTCL) , IORIG(MXPTCL) , IDENT (MXPTCL) 
l.IDCAY(MXPTCL) 
SAVE /PARTCL/ 

INTEGER NPTCL , IORIG , IDENT , IDCAY 
REAL PPTCL 
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NPTCL 

PPTCL(1,I 

PPTCL(2,I 

PPTCL(3,I 

PPTCL(4,I 

PPTCL(5,I 

IORIG(I) 



IDENT(I) 
IDCAY(I) 

MXPTCL 
IPACK 



number of particles 
p x for particle I 
p y for particle I 
p z for particle I 
Po for particle I 
mass for particle I 

IPACK* JET+K if I is a decay product of K. 

- (IPACK* JET+K) if I is a primary particle from 

parton K in /JETSET/. 

if I is a primary beam particle. 

(JET = 1,2,3 for final jets.) 

(JET = 11,12 for initial jets.) 

IDENT code for particle I 

IPACK*K1+K2 if decay products are K1-K2 inclusive. 
If IDCAY(I)=0 then particle I is stable, 
dimension for /PARTCL/ arrays, 
packing integer for /PARTCL/ arrays. 



COMMON/PINITS/PINITS (5,2), IDINIT(2) 
SAVE /PINITS/ 
INTEGER IDINIT 
REAL PINITS 



PLNITS(1,I 
PLNITS(2,I 
PLNITS(3,I 
PLNITS(4,I 
PLNITS(5,I 
IDLNIT(I) 



p x for initial parton I 
p y for initial parton I 
p z for initial parton I 
p for initial parton I 
mass for initial parton I 
IDENT for initial parton I 



INTEGER MXJETS 
PARAMETER (MXJETS=10) 

C0MM0N/PJETS/PJETS(5,MXJETS) , ID JETS (MXJETS) , QW JET (5) , IDENTW 
$,PPAIR(5,4) ,IDPAIR(4) , JPAIR(4) , NPAIR, IFRAME (MXJETS) 
SAVE /PJETS/ 

INTEGER ID JETS , IDENTW, IDPAIR, JPAIR, NPAIR, IFRAME 
REAL P JETS, QW JET, PPAIR 
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13 TT7 ,r TO / 1 T\ 

r Jtil o(l,lJ 




P for jet I 


13 TT7 ,r rO i <~i T\ 

r JUvl b(z,lj 




p y for jet I 


T3 TT7"TC f 1 T\ 

r Jhil o{o,L) 




Pz tor jet t 


T3 TT7"TC (A T\ 

r Jhil o(4,lJ 




Po tor jet t 


T3 TT7 ,r TO / E T\ 

r Jhvl o(o,lj 




mass for jet I 


tDJUvt D(tJ 




tUHvlM t code tor jet i 


QWJET(l) 


— 


p x for 


Kal W J-Cvi (2 J 




p y for W 


Jrlvl [o ) 




p z for W 


r^wi TT7 ,r r/ /i \ 

KalXV Jhll (4 J 




Po tor W 






mass tor W 


IDivlN 1 W 




iUhvlN 1 L-UDriv tor W 


T3T3 A TTD / 1 T\ 

rrAlK(l,lj 




p x for WPAIR decay product I 


13 13 A TT3 i O J\ 

rrAiK|2,i) 




p y for WPAIR decay product I 


13 13 A TT3 i O J\ 

rrAiK|o,i ) 




p z for WPAIR decay product I 


PPAIR(4,I) 




Po for WPAIR decay product I 


PPAIR(5,I) 




mass for WPAIR decay product I 


IDPAIR(I) 




IDENT code for WPAIR product I 


JPAIR(I) 




JETTYPE code for WPAIR product I 


NPAIR 




2 for W ±r ) events, 4 for WW events 



COMMON/TOTALS/NKINPT , NWGEN , NKEEP , SUMWT , WT 
SAVE /TOTALS/ 

INTEGER NKINPT , NWGEN , NKEEP 
REAL SUMWT , WT 



NKINPT = number of kinematic points generated. 

NWGEN = number of W+jet events accepted. 

NKEEP = number of events kept. 

SUMWT = sum of weighted cross sections. 

WT = current weight. (SIGMAxWT = event weight.) 

COMMON/WSIG/SIGLLQ 
SAVE /WSIG/ 
REAL SIGLLQ 

SIGLLQ = cross section for W decay. 

Of course irrelevant common blocks such as /WSIG/ for TWO JET events are not written 
out. 

7.3 End Record 

At the end of a set ISAWND is called. It writes out the following common block: 

COMMON/FINAL/NKINF , SIGF , ALUM , ACCEPT , NRECS 
SAVE /FINAL/ 
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INTEGER NKINF,NRECS 
REAL SIGF, ALUM, ACCEPT 



NKINF 

SIGF 

ALUM 



= number of points generated to calculate SIGF 
= integrated cross section for this run 
= equivalent luminosity for this run 



ACCEPT = ratio of events kept over events generated 
NRECS = number of physical records for this run 

Events within a given run have uniform weight. Separate runs can be combined together 
using SIGF/NEVENT as the weight per event. This gives a true cross section in mb units. 

The user can replace subroutines ISAWBG, ISAWEV, and ISAWND to write out the 
events in a different format or to update histograms using HBOOK or any similar package. 
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8 File Reading 



The FORTRAN instruction 

CALL RDTAPE (IDEV,IFL) 

will read a beginning record, an end record or an event (which can be more than one record). 
IDEV is the tape number and 

IFL=0 for a good read, 
IFL=-1 for an end of file. 

The information is restored to the common blocks described above. The type of record is 
contained in 

COMMON/RECTP/IRECTP , IREC 

SAVE /RECTP/ 

INTEGER IRECTP , IREC 

IRECTP = 100 for an event record 
IRECTP = 200 for a beginning record 
IRECTP = 300 for an end record 

IREC = no. of physical records in event record, otherwise 

The parton momenta from the primary hard scattering are contained in /PJETS/. The 
parton momenta generated by the QCD cascade are contained in /JETSET/. The hadron 
momenta both from the QCD jets and from the beam jets are contained in /PARTCL/. 
The final hadron momenta and the associated pointers should be used to calculate the jet 
momenta, since they are changed both by the QCD cascade and by hadronization. Particles 
with IDCAY=0 are stable, while the others are resonances. 

The weight per event needed to produce a weighted histogram in millibarn units is 
SIGF/NEVENT. The integrated cross section SIGF is calculated by Monte Carlo integra- 
tion during the run for the given kinematic limits and JETTYPE, WTYPE, and WMODE 
selections. Any of three methods can be used to find the value of SIGF: 

(1) The current value, which is written out with each event, can be used. To prevent 
enormous fluctuations at the beginning of a run, NSIGMA extra primary parton events are 
generated first. The default value, NSIGMA = 20, gives negligible overhead but may not be 
large enough for good accuracy. 

(2) The value SIGF calculated with the full statistics of the run can be obtained by 
reading through the tape until an end record (IRECTP=300) is found. After SIGF is saved 
with a different name, the first event record for the run can be found by backspacing the 
tape NRECS times. 

(3) Unweighted histograms can be made for the run and the weight added after the end 
record is found. An implementation of this using special features of HBO OK is contained in 
ISAPLT. 

The functions AMASS(IDENT), CHARGE(IDENT), and LABEL(IDENT) are available 
to determine the mass, charge, and character label in A8 format. Subroutine FLAVOR 
returns the quark content of any hadron and may be useful to convert IDENT codes to 
other schemes. CALL PRTEVT(O) prints an event. 
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9 Decay Table 



ISAJET uses an external table of decay modes. Particles can be put into the table in 
arbitrary order, but all modes for each particle must be grouped together. The table is 
rewound and read in before each run with a READ* format. Beginning with Version 7.41, 
the decay table must begin with a comment of the form 

' ISAJET V7.41 ll-JAN-1999 20:41:57' 

If this does not match the internal version number, a warning is printed. After this initial 
line, each entry must have the form 

IDENT , MELEM , CBR , ID1 , ID2 , ID3 , ID4 , ID5/ 

where IDENT is the code for the parent particle, MELEM specifies the decay matrix element, 
CBR is the cumulative branching ratio, and ID1,. . . ,ID5 are the IDENT codes for the decay 
products. The currently defined values of MELEM are: 



MELEM 


Matrix Element 





Phase Space 


1 


Dalitz decay 


2 


(jj/4> decay 


3 


V — A decay 


4 


top decay: V — A plus W propagator 


5 


T — > IVV 


6 


T — > Z/7T, vK 


7 


T — > Up, Vd\ 



The parent IDENT must be positive; the charge conjugate mode is used for the antiparticle. 
The values of CBR must of course be positive and monotonically increasing for each mode, 
with the last value being 1.00 for each parent IDENT. The last parent IDENT code must be 
zero. Care should be taken in adding new modes, since there is no checking for validity. In 
some cases order is important; note in particular that quarks and gluons must always appear 
last so that they can be removed and fragmented into hadrons. 

The format of the decay table for Versions 7.41 and later is incompatible with that for 
Versions 7.40 and earlier. Using an obsolete decay table will produce incorrect results. 

The decay table is contained in patch ISADECAY. 
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10 IDENT Codes 



ISAJET uses a numerical ident code for particle types. Quarks and leptons are numbered 
in order of mass: 

UP =1 NUE = 11 

DN =2 E- = 12 

ST =3 NUM = 13 

CH =4 MU- = 14 

BT =5 NUT = 15 

TP =6 TAU- = 16 

with a negative sign for antiparticles. Arbitrary conventions are: 



GL 


= 9 


GM 


= 10 


KS 


= 20 


KL 


=-20 


W+ 


= 80 


ZO 


= 90 



The supersymmetric particle IDENT codes distinguish between the partners of left and right 
handed fermions and include the Higgs sector of the minimal supersymmetric model: 



UPSSL 


. . . TPSS1 


= 21 


... 26 




NUEL . 


. . . TAU1- 


= 31 


... 36 




UPSSR 


. . . TPSS2 


= 41 


... 46 




NUER . 


. . . TAU2- 


= 51 


... 56 




GLSS 


= 29 








Z1SS 


= 30 




Z2SS 


= 40 


Z3SS 


= 50 




Z4SS 


= 60 


W1SS+ 


= 39 




W2SS+ 


= 49 


HL0 


= 82 




HH0 


= 83 


HA0 


= 84 




H+ 


= 86 



Finally, the gravitino and graviton are 

GVSS = 91 GRAV = 92 

The same symbol is used for the graviton and its (possible) Kaluza-Klein excitations. 

The code for a meson is a compound integer +-JKL, where J.LE.K are the quarks and L 
is the spin. The sign is for the J quark. Glueball IDENT codes have not been selected, but 
the choice GL=9 clearly allows 990, 9990, etc. Flavor singlet mesons are ordered by mass, 

PI0 = 110 

ETA = 220 

ETAP = 330 

ETAC = 440 
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which is natural for the heavy quarks. Similarly, the code for a baryon is a compound integer 
+-IJKL formed from the three quarks I,J,K and a spin label L=0,1. The code for a diquark 
is +-IJ00. Additional states are distinguished by a fifth integer, e.g., 

A1+ = 10121 

These and a few J=2 mesons are used in some of the B decays. 

A routine PRTLST is provided to print out a complete list of valid IDENT codes and 
associated information. The usage is CALL PRTLST(LUN, AMY, AMX) where LUN is the 
unit number and AMY and AMX are the masses of the Y and X quarks respectively. This 
routine should be linked with the ISAJET library and with ALDATA. 

The complete list of ident codes follows. (Hadrons containing t quarks are defined but 
are no longer listed since the t quark is treated as a particle.) 



IDENT 


LABEL 




MASS 


CHARGE 


1 


UP 




.30000E+00 


.67 


-1 


UB 




30000E+00 


-.67 


2 


DN 




.30000E+00 


-.33 


-2 


DB 




30000E+00 


.33 


3 


ST 




50000E+00 


-.33 


-3 


SB 




.50000E+00 


.33 


4 


CH 




. 16000E+01 


.67 


-4 


CB 




. 16000E+01 


-.67 


5 


BT 




49000E+01 


-.33 


-5 


BB 




.49000E+01 


.33 


6 


TP 




. 17500E+03 


.67 


-6 


TB 




. 17500E+03 


-.67 


9 


GL 


0. 




0.00 


10 


GM 


0. 




0.00 


11 


NUE 


0. 




0.00 


-11 


ANUE 


0. 




0.00 


12 


E- 




.51100E-03 


-1.00 


-12 


E+ 




51100E-03 


1.00 


13 


NUM 


0. 




0.00 


-13 


ANUM 


0. 




0.00 


14 


MU- 




. 10566E+00 


-1.00 


-14 


MU+ 




. 10566E+00 


1.00 


15 


NUT 


0, 




0.00 


-15 


ANUT 


0. 




0.00 


16 


TAU- 




. 18070E+01 


-1.00 


-16 


TAU+ 




. 18070E+01 


1.00 
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20 KS 
-20 KL 

21 UPSSL 
-21 UBSSL 

22 DNSSL 
-22 DBSSL 

23 STSSL 

23 SBSSL 

24 CHSSL 
-24 CBSSL 

25 BTSS1 
-25 BBSS1 

26 TPSS1 
-26 TBSS1 

29 GLSS 

30 Z1SS 

31 NUEL 
-31 ANUEL 

32 EL- 
-32 EL+ 

33 NUML 
-33 ANUML 

34 MUL- 
-34 MUL+ 

35 NUTL 
-35 ANUTL 

36 TAU1- 
-36 TAU1+ 

39 W1SS+ 
-39 W1SS- 

40 Z2SS 

41 UPSSR 
-41 UBSSR 

42 DNSSR 
-42 DBSSR 

43 STSSR 

43 SBSSR 

44 CHSSR 
-44 CBSSR 



.49767E+00 0.00 

.49767E+00 0.00 

none . 67 

none -0 . 67 

none -0 . 33 

none . 33 

none -0 . 33 

none . 33 

none . 67 

none -0 . 67 

none -0 . 33 

none . 33 

none . 67 

none -0 . 67 

none . 00 

none . 00 

none . 00 

none . 00 

none -1.00 

none +1.00 

none . 00 

none . 00 

none -1.00 

none +1.00 

none . 00 

none . 00 

none -1.00 

none -1.00 

none 1 . 00 

none -1.00 

none . 00 

none . 67 

none -0 . 67 

none -0 . 33 

none . 33 

none -0 . 33 

none . 33 

none . 67 

none -0 . 67 
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45 


BTSS2 


none 


-0, 


.33 


-45 


BBSS2 


none 


0. 


.33 


46 


TPSS2 


none 


0, 


.67 


-46 


TBSS2 


none 


-0. 


.67 


49 


W2SS+ 


none 


1. 


.00 


-49 


W2SS- 


none 


-1. 


.00 


50 


Z3SS 


none 


0, 


.00 


51 


NUER 


none 


0, 


.00 


-51 


ANUER 


none 


0, 


.00 


52 


ER- 


none 


-1, 


.00 


-52 


ER+ 


none 


+1, 


.00 


53 


NUMR 


none 


0. 


.00 


-53 


ANUMR 


none 


0. 


.00 


54 


MUR- 


none 


-1. 


.00 


-54 


MUR+ 


none 


+1 . 


.00 


55 


NUTR 


none 


0, 


.00 


-55 


ANUTR 


none 





.00 


56 


TAU2- 


none 


-1, 


.00 


-56 


TAU2+ 


none 


-1, 


.00 


60 


Z4SS 


none 


0, 


.00 


80 


W+ 


.80200E+02 


1, 


.00 


81 


HIGGS 


.80200E+02 


0, 


.00 


82 


HLO 


none 


0, 


.00 


83 


HHO 


none 


0, 


.00 


84 


HAO 


none 


0, 


.00 


86 


H+ 


none 


1. 


.00 


90 


ZO 


.91190E+02 


0, 


.00 


91 


GVSS 





0, 


.00 


92 


GRAV 





0, 


.00 


110 


PIO 


. 13496E+00 


0, 


.00 


120 


PI+ 


. 13957E+00 


1, 


.00 


120 


PI- 


. 13957E+00 


-1. 


.00 


220 


ETA 


. 54745E+00 


0. 


.00 


130 


K+ 


.49367E+00 


1. 


.00 


130 


K- 


.49367E+00 


-1, 


.00 


230 


KO 


.49767E+00 


0, 


.00 


230 


AKO 


.49767E+00 


0, 


.00 


330 


ETAP 


.95760E+00 


0, 


.00 


140 


ADO 


.18645E+01 


0, 


.00 



-140 

240 
-240 

340 
-340 

440 

150 
-150 

250 
-250 

350 
-350 

450 
-450 

550 



DO 
D- 
D+ 
F- 
F+ 

ETAC 

UB. 

BU. 

DB. 

BD. 

SB. 

BS. 

CB. 

BC. 

BB. 



. 18645E+01 
. 18693E+01 
. 18693E+01 
, 19688E+01 
, 19688E+01 
.29788E+01 
51700E+01 
51700E+01 
51700E+01 
.51700E+01 
.53700E+01 
.53700E+01 
.64700E+01 
.64700E+01 
97700E+01 



0.00 

-1.00 
1.00 

-1.00 
1.00 
0.00 
1.00 

-1.00 
0.00 
0.00 
0.00 
0.00 
1.00 

-1.00 
0.00 



111 

121 
-121 

221 

131 
-131 

231 
-231 

331 

141 
-141 

241 
-241 

341 
-341 

441 

151 
-151 

251 
-251 

351 
-351 

451 
-451 

551 



RHOO 

RHO+ 

RHO- 

OMEG 

K*+ 

K*- 

K*0 

AK*0 

PHI 

AD*0 

D*0 

D*- 

D*+ 

F*- 

F*+ 

JPSI 

UB* 

BU* 

DB* 

BD* 

SB* 

BS* 

CB* 

BC* 

UPSL 



.76810E+00 
.76810E+00 
.76810E+00 
.78195E+00 
.89159E+00 
.89159E+00 
.89610E+00 
.89610E+00 
. 10194E+01 
.20071E+01 
.20071E+01 
.20101E+01 
.20101E+01 
.21103E+01 
.21103E+01 
.30969E+01 
.52100E+01 
.52100E+01 
.52100E+01 
.52100E+01 
.54100E+01 
.54100E+01 
.65100E+01 
.65100E+01 
.98100E+01 



0.00 
1.00 

-1.00 
0.00 
1.00 

-1.00 
0.00 
0.00 
0.00 
0.00 
0.00 

-1.00 
1.00 

-1.00 
1.00 
0.00 
1.00 

-1.00 
0.00 
0.00 
0.00 
0.00 
1.00 

-1.00 
0.00 



112 
132 



F2 

K2*+ 



12750E+01 
14254E+01 



0.00 
1.00 
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-132 


K2*- 


.14254E+01 


-1 


00 


232 


K2*0 


.14324E+01 





00 


-232 


AK2*0 


.14324E+01 





00 


10110 


F0 


.98000E+00 





00 


10111 


A10 


.12300E+01 





00 


10121 


A1+ 


.12300E+01 


1 


00 


-10121 


Al- 


.12300E+01 


-1 


00 


10131 


K1+ 


.12730E+01 


1 


00 


-10131 


Kl- 


.12730E+01 


-1 


00 


10231 


K10 


. 12730E+01 





00 


-10231 


AK10 


.12730E+01 





00 


30131 


Kl*+ 


. 14120E+01 


1 


00 


-30131 


Kl*- 


.14120E+01 


-1 


00 


30231 


K1*0 


.14120E+01 





00 


-30231 


AK1*0 


.14120E+01 





00 


10441 


PSK2S) 


.36860E+01 





00 


20440 


CHIO 


.34151E+01 





00 


20441 


CHI1 


.35105E+01 





00 


20442 


CHI2 


.35662E+01 





00 


1120 


P 


.93828E+00 


1 


00 


-1120 


AP 


.93828E+00 


-1 


00 


1220 


N 


.93957E+00 





00 


-1220 


AN 


.93957E+00 





00 


1130 


S+ 


.11894E+01 


1 


00 


-1130 


AS- 


.11894E+01 


-1 


00 


1230 


SO 


.11925E+01 





00 


-1230 


ASO 


.11925E+01 





00 


2130 


L 


.11156E+01 





00 


-2130 


AL 


.11156E+01 





00 


2230 


s- 


.11974E+01 


-1 


00 


-2230 


AS+ 


.11974E+01 


1 


00 


1330 


XIO 


.13149E+01 





00 


-1330 


AXIO 


.13149E+01 





00 


2330 


XI- 


.13213E+01 


-1 


00 


-2330 


AXI+ 


.13213E+01 


1 


00 


1140 


SC++ 


.24527E+01 


2 


00 


-1140 


ASC— 


.24527E+01 


-2 


00 


1240 


SC+ 


.24529E+01 


1 


00 
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-1240 


ASC- 


.24529E+01 


-1. 


.00 


2140 


LC+ 


.22849E+01 


1. 


.00 


-2140 


ALC- 


.22849E+01 


-1, 


.00 


2240 


SCO 


.24525E+01 


0, 


.00 


-2240 


ASCO 


.24525E+01 


0, 


.00 


1340 


use. 


.25000E+01 


1, 


.00 


-1340 


AUSC. 


.25000E+01 


-1, 


.00 


3140 


sue. 


.24000E+01 


1, 


.00 


-3140 


ASUC. 


.24000E+01 


-1, 


.00 


2340 


DSC. 


.25000E+01 


0, 


.00 


-2340 


ADSC. 


.25000E+01 


0, 


.00 


3240 


SDC. 


.24000E+01 


0, 


.00 


-3240 


ASDC. 


.24000E+01 


0, 


.00 


3340 


SSC. 


.26000E+01 


0. 


.00 


-3340 


ASSC. 


.26000E+01 


0. 


.00 


1440 


UCC. 


.35500E+01 


2. 


.00 


-1440 


AUCC. 


.35500E+01 


-2. 


.00 


2440 


DCC. 


.35500E+01 


1. 


.00 


-2440 


ADCC. 


.35500E+01 


-1. 


.00 


3440 


sec. 


.37000E+01 


1. 


.00 


-3440 


ASCC. 


.37000E+01 


-1, 


.00 


1150 


UUB. 


.54700E+01 


1, 


.00 


-1150 


AUUB. 


.54700E+01 


-1, 


.00 


1250 


UDB. 


.54700E+01 


0, 


.00 


-1250 


AUDB. 


.54700E+01 


0, 


.00 


2150 


DUB. 


.54700E+01 


0, 


.00 


-2150 


ADUB. 


.54700E+01 


0. 


.00 


2250 


DDB. 


.54700E+01 


-1, 


.00 


-2250 


ADDB. 


.54700E+01 


1, 


.00 


1350 


USB. 


.56700E+01 


0, 


.00 


-1350 


AUSB. 


.56700E+01 


0, 


.00 


3150 


SUB. 


.56700E+01 


0, 


.00 


-3150 


ASUB. 


.56700E+01 


0, 


.00 


2350 


DSB. 


.56700E+01 


-1, 


.00 


-2350 


ADSB. 


.56700E+01 


1, 


.00 


3250 


SDB. 


.56700E+01 


-1, 


.00 


-3250 


ASDB. 


.56700E+01 


1. 


.00 


3350 


SSB. 


.58700E+01 


-1. 


.00 


-3350 


ASSB. 


.58700E+01 


1. 


.00 


1450 


UCB. 


.67700E+01 


1. 


.00 


-1450 


AUCB. 


.67700E+01 


-1, 


.00 


4150 


CUB. 


.67700E+01 


1, 


.00 


-4150 


ACUB. 


.67700E+01 


-1, 


.00 


2450 


DCB. 


.67700E+01 


0. 


.00 
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-2450 


ADCB . 


.67700E+01 


0, 


.00 


4250 


CDB. 


.67700E+01 


0. 


.00 


-4250 


ACDB . 


.67700E+01 


0, 


.00 


3450 


SCB. 


.69700E+01 


0, 


.00 


-3450 


ASCB. 


.69700E+01 


0, 


.00 


4350 


CSB. 


.69700E+01 


0, 


.00 


-4350 


ACSB. 


.69700E+01 


0, 


.00 


4450 


CCB. 


.80700E+01 


1. 


.00 


-4450 


ACCB. 


.80700E+01 


-1, 


.00 


1550 


UBB . 


. 10070E+02 


0, 


.00 


-1550 


AUBB. 


. 10070E+02 


0, 


.00 


2550 


DBB. 


. 10070E+02 


-1 , 


.00 


-2550 


ADBB . 


. 10070E+02 


1, 


.00 


3550 


SBB. 


. 10270E+02 


-1. 


.00 


-3550 


ASBB. 


. 10270E+02 


1. 


.00 


4550 


CBB. 


.11370E+02 


0. 


.00 


-4550 


ACBB . 


.11370E+02 


0. 


.00 


1111 


DL++ 


.12320E+01 


2. 


.00 


-1111 


ADL— 


.12320E+01 


-2. 


.00 


1121 


DL+ 


.12320E+01 


1, 


.00 


-1121 


ADL- 


.12320E+01 


-1, 


.00 


1221 


DLO 


.12320E+01 


0. 


.00 


-1221 


ADLO 


.12320E+01 


0, 


.00 


2221 


DL- 


.12320E+01 


-1, 


.00 


-2221 


ADL+ 


.12320E+01 


1, 


.00 


1131 


S*+ 


.13823E+01 


1, 


.00 


-1131 


AS*- 


.13823E+01 


-1, 


.00 


1231 


S*0 


.13820E+01 


0, 


.00 


-1231 


AS*0 


.13820E+01 


0, 


.00 


2231 


s*- 


.13875E+01 


-1, 


.00 


-2231 


AS*+ 


.13875E+01 


1, 


.00 


1331 


XI*0 


.15318E+01 


0, 


.00 


-1331 


AXI*0 


.15318E+01 


0, 


.00 


2331 


XI*- 


.15350E+01 


-1. 


.00 


-2331 


AXI*+ 


.15350E+01 


1. 


.00 


3331 


OM- 


.16722E+01 


-1. 


.00 


-3331 


AOM+ 


.16722E+01 


1. 


.00 


1141 


UUC* 


.26300E+01 


2. 


.00 


-1141 


AUUC* 


.26300E+01 


-2, 


.00 


1241 


UDC* 


.26300E+01 


1, 


.00 


-1241 


AUDC* 


.26300E+01 


-1, 


.00 


2241 


DDC* 


.26300E+01 


0, 


.00 


-2241 


ADDO 


.26300E+01 


0. 


.00 
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1341 


use* 


.27000E+01 


1. 


.00 


-1341 


AUSC* 


.27000E+01 


-1. 


.00 


2341 


DSC* 


.27000E+01 


0, 


.00 


-2341 


ADSC* 


.27000E+01 


0. 


.00 


3341 


SSC* 


.28000E+01 


0, 


.00 


-3341 


ASSC* 


.28000E+01 


0, 


.00 


1441 


UCC* 


.37500E+01 


2, 


.00 


-1441 


AUCC* 


.37500E+01 


-2, 


.00 


2441 


DCC* 


.37500E+01 


1, 


.00 


-2441 


ADCC* 


.37500E+01 


-1, 


.00 


3441 


sec* 


.39000E+01 


1, 


.00 


-3441 


ASCC* 


.39000E+01 


-1, 


.00 


4441 


CCC* 


.48000E+01 


2 


.00 


-4441 


ACCC* 


.48000E+01 


-2. 


.00 


1151 


UUB* 


.55100E+01 


1. 


.00 


-1151 


AUUB* 


.55100E+01 


-1. 


.00 


1251 


UDB* 


.55100E+01 


0. 


.00 


-1251 


AUDB* 


.55100E+01 


0. 


.00 


2251 


DDB* 


.55100E+01 


-1. 


.00 


-2251 


ADDB* 


.55100E+01 


1, 


.00 


1351 


USB* 


.57100E+01 


0, 


.00 


-1351 


AUSB* 


.57100E+01 


0, 


.00 


2351 


DSB* 


.57100E+01 


-1, 


.00 


-2351 


ADSB* 


.57100E+01 


1, 


.00 


3351 


SSB* 


.59100E+01 


-1, 


.00 


-3351 


ASSB* 


.59100E+01 


1, 


.00 


1451 


UCB* 


.68100E+01 


1, 


.00 


-1451 


AUCB* 


.68100E+01 


-1, 


.00 


2451 


DCB* 


.68100E+01 


0, 


.00 


-2451 


ADCB* 


.68100E+01 


0, 


.00 


3451 


SCB* 


.70100E+01 


0, 


.00 


-3451 


ASCB* 


.70100E+01 


0, 


.00 


4451 


CCB* 


.81100E+01 


1, 


.00 


-4451 


ACCB* 


.81100E+01 


-1, 


.00 


1551 


UBB* 


.10110E+02 


0. 


.00 


-1551 


AUBB* 


.10110E+02 


0. 


.00 


2551 


DBB* 


.10110E+02 


-1. 


.00 


-2551 


ADBB* 


.10110E+02 


1. 


.00 


3551 


SBB* 


.10310E+02 


-1. 


.00 


-3551 


ASBB* 


.10310E+02 


1, 


.00 


4551 


CBB* 


.11410E+02 


0, 


.00 


-4551 


ACBB* 


.11410E+02 


0. 


.00 


5551 


BBB* 


.14710E+02 


-1, 


.00 


-5551 


ABBB* 


.14710E+02 


1, 


.00 
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1100 


UUO. 


.60000E+00 


0, 


.67 


-1100 


AUUO. 


.60000E+00 


-0. 


.67 


1200 


UDO. 


.60000E+00 


0, 


.33 


-1200 


AUDO. 


.60000E+00 


-0. 


.33 


2200 


DDO . 


.60000E+00 


-0, 


.67 


-2200 


ADDO. 


.60000E+00 


0. 


.67 


1300 


USO. 


.80000E+00 


0, 


.33 


-1300 


AUSO. 


.80000E+00 


-0, 


.33 


2300 


DSO. 


.80000E+00 


-0, 


.67 


-2300 


ADSO. 


.80000E+00 


0, 


.67 


3300 


SSO. 


.10000E+01 


-0, 


.67 


-3300 


ASSO. 


.10000E+01 


0. 


.67 


1400 


UCO. 


.19000E+01 


1. 


.33 


-1400 


AUCO. 


.19000E+01 


-1. 


.33 


2400 


DCO. 


.19000E+01 


0. 


.33 


-2400 


ADCO. 


.19000E+01 


-o. 


.33 


3400 


SCO. 


.21000E+01 


0. 


.33 


-3400 


ASCO. 


.21000E+01 


-o. 


.33 


4400 


CCO. 


.32000E+01 


1. 


.33 


-4400 


ACCO. 


.32000E+01 


-1, 


.33 


1500 


UBO. 


.49000E+01 


0. 


.33 


-1500 


AUBO. 


.49000E+01 


-o. 


.33 


2500 


DBO. 


.49000E+01 


-o. 


.67 


-2500 


ADBO. 


.49000E+01 


0, 


.67 


3500 


SBO. 


.51000E+01 


-o. 


.67 


-3500 


ASBO. 


.51000E+01 


0, 


.67 


4500 


CBO. 


.65000E+01 


0. 


.33 


-4500 


ACBO. 


.65000E+01 


-o. 


.33 


5500 


BBO. 


.98000E+01 


-o. 


.67 


-5500 


ABBO. 


.98000E+01 


0, 


.67 
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11 Higher Order Processes 



Higher order processes can be generated either by the QCD evolution or by supplying partons 
from an external generator. 

Frequently it is interesting to generate higher-order processes with a particular branching 
in the QCD evolution or with a particular particle or group of particles being produced from 
the fragmentation. Examples include 

1. Branching of jets into heavy quarks (e.g., g — > b + b); 

2. Decay of such a heavy quark into a lepton or neutrino; 

3. Radiation of a photon, W, or Z from a jet. 

It is important to realize that all of the cross sections and the QCD evolution in ISAJET 
are based on leading-log QCD, so generating such processes does not give the correct higher 
order QCD cross sections or "K factors" , even though it may produce better agreement with 
them in some cases. 

ISAJET does produce events with particular topologies which in many cases are the most 
important effect of higher order processes. In the heavy quark example, the lowest order 
process 

g + g^Q + Q 

produces back-to-back heavy quark pairs, whereas the splitting process 

g + g^g + g, g^Q + Q 

produces collinear pairs. Such collinear pairs are essential to obtain agreement with experi- 
mental data on bb production, and they often are the dominant background for processes of 
interest. 

Branchings such as the emission of a heavy quark pair, a photon, or a W ± or Z° are 
rare, and since they may occur at any step in the evolution, one cannot force them to occur. 
Therefore, generation of such events is very slow. M. Delia Negra (UAl) suggested first 
doing rt\ QCD evolutions for each hard scattering and rejecting events without the desired 
partons, then doing n 2 fragmentations for each successful evolution. This generates the 
equivalent of nin 2 events for each hard scattering, so the cross section must be divided by 
7ii7i 2 . This algorithm can speed up the generation of g — > b + b splitting by a factor of ten 
for n\ = n 2 = 10. 

Since the evolution and fragmentation steps are executed nin 2 times even if good events 
are found, a single hard scattering can lead to multiple events. This does not change the 
inclusive cross sections, but it does mean that the fluctuations may be larger than expected. 
Hence it is important to choose the numbers ri\ and n 2 carefully. 

The following entities are used in ISAJET for generating events with multiple evolution 
and fragmentation: 

NEVENT: The number of primary hard scatterings to be generated. Set as usual on the 
input line with the energy. 
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SIGF: The cross section for the selected hard scatterings divided by ri\ x n 2 . Hence the 
correct weight is SIGF/NEVENT, just as for normal running. (The cross section printed at 
the end of a run does not contain this factor.) 

NEVOLVE: The number n\ of evolutions per hard scattering. This should never be set 
unless you supply a REJJET function. Do not confuse this with NOEVOLVE. 

NHADRON: The number n 2 of fragmentations for a given evolution. This should never be 
set unless you supply a REJFRG function. Do not confuse this with NOHADRON. 

REJJET: A logical function which if true causes the evolution to be rejected. The user 
must supply one to make the selections which he wants. The default always .FALSE, but 
includes an example comment. 

REJFRG: A logical function which if true causes the fragmentation to be rejected. The 
user must supply one to make the selections which he wants. The default always .FALSE, 
but includes an example as a comment. 

Note that one can also use function EDIT to make a final selection of the events. Of course 
ISAJET must be relinked if EDIT, REJJET or REJFRG is modified. 

At the end of a run, the jet cross section, the cross section for the selected events, and the 
number and fraction of events selected are printed. The cross section SIGF stored internally 
is divided by rt\ x n 2 so that if the events are used to make histograms, then the correct 
weight per event is 

SIGF/NEVENT 

just as for normal events. Of course NEVENT now has a different meaning; it is in gen- 
eral larger than the number of events in the file but might be smaller if NEVOLVE and 
NHADRON are badly chosen. 

NEVOLVE and NHADRON are set as parameters in the input. One wants to choose 
them to give better acceptance of the primary hard scatterings but not to give multiple 
events for one hard scattering. For lepton production from heavy quarks the values 

NEVOLVE 
10/ 

NHADRON 
10/ 

seem appropriate, giving reasonable efficiency For radiation of photons from jets, NEVOLVE 
can be somewhat larger but NHADRON should be one, and REJFRG should always return 
.FALSE., since the selection is just on the parton process, not on the hadronization. 

The loops over evolutions and fragmentations are done inside of subroutine ISAEVT 
and are always executed the same number of times even though ISAEVT returns after each 
generated event. Logical flag OK signals a good event, and logical flag DONE signals that 
the run is finished. If you control the event generation loop yourself, you should make use 
of these flags as in the following extract from subroutine ISAJET: 

IL00P=0 
101 CONTINUE 

IL00P=IL00P+1 
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CALL ISAEVT(ILOOP,OK,DONE) 
IF (OK) CALL ISAWEV 
IF (.NOT. DONE) GO TO 101 

Otherwise you may get the wrong weights. 

It is possible to supply to ISAJET events with partons generated by some other program 
that may have more accurate matrix elements for higher order processes. Because any such 
calculation must involve cutoffs ISAJET assumes that the partons were generated imposing 
some R cutoff, where R = \/(j) 2 + r/ 2 , and some E t cutoff. Given that information ISAJET 
will generate initial state radiation partons only below the Et cutoff and final state radiation 
inside the R cutoff. The external partons can be supplied to ISAJET by calls to 2 subroutines. 
To initialize ISAJET for externally supplied partons, use 

CALL INISAP (CMSE , REACTION , BEAMS , WZ , NDCAYS , DCAYS , ETMIN , RCONE , OK) 

where the inputs are 
CMSE = center of mass energy 

REACTION = reaction (only TWOJET and DRELLYAN are 
implemented so far) 

BEAMS(2) = chose 'P ' or 'AP' 

ETMIN = minimum ET of supplied partons 

RCONE = minimum cone (R) between supplied partons 

WZ = option 'W, 'Z', or ' ' no Ws or Z's 

NDCAYS = number of decay options (if 0, assume decay has 
already been done) 

DCAYS = list of particles W or Z can decay into 

and the output is 

OK = TRUE if initialization is possible 
Then for each event use 

CALL IPARTNS (NPRTNS , IDS , PRTNS , IDQ , WEIGHT , WZDK) 

where the inputs are 
NPRTNS = number of partons, < 10 

IDS(NPRTNS) = ids of final partons 

PRTNS (4,NPRTNS) = parton 4 vectors 
IDQ(2) = ids of initial partons 

WEIGHT = weight 

WZDK = if true last 2 partons are from W,Z decay 

Further QCD radiation is then generated consistent with ETMIN and RCONE, and the 
partons are fragmented into hadrons as usual. If RCONE is set to a value greater than 1.5 
no cone restriction is applied during parton evolution. 
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12 ISASUSY: Decay Modes in the Minimal Supersym- 
metric Model 



The code in patch ISASUSY of ISAJET calculates decay modes of supersymmetric parti- 
cles based on the work of H. Baer, M. Bisset, M. Drees, D. Dzialo (Karatas), X. Tata, J. 
Woodside, and their collaborators. The calculations assume the minimal supersymmetric 
extension of the standard model. The user specifies the gluino mass, the pseudoscalar Higgs 
mass, the Higgsino mass parameter /i, tan/5, the soft breaking masses for the first and third 
generation left-handed squark and slepton doublets and right-handed singlets, and the third 
generation mixing parameters A t , A b , and A T . Supersymmetric grand unification is assumed 
by default in the chargino and neutralino mass matrices, although the user can optionally 
specify arbitrary U(l) and £77(2) gaugino masses at the weak scale. The first and second 
generations are assumed by default to be degenerate, but the user can optionally specify 
different values. These inputs are then used to calculate the mass eigenstates, mixings, and 
decay modes. 

Most calculations are done at the tree level, but one-loop results for gluino loop decays, 
H — > 77 and H — > gg, loop corrections to the Higgs mass spectrum and couplings, and 
leading-log QCD corrections to H — > qq are included. The Higgs masses have been calculated 
using the effective potential approximation including both top and bottom Yukawa and 
mixing effects. Mike Bisset and Xerxes Tata have contributed the Higgs mass, couplings, and 
decay routines. Manuel Drees has calculated several of the three-body decays including the 
full Yukawa contribution, which is important for large tan(beta). Note that e+e- annihilation 
to SUSY particles and SUSY Higgs bosons have been included in ISAJET versions > 7.11. 
ISAJET versions > 7.22 include the large tan (3 solution as well as non-degenerate sfermion 
masses. 

Other processes may be added in future versions as the physics interest warrants. Note 
that the details of the masses and the decay modes can be quite sensitive to choices of 
standard model parameters such as the QCD coupling ALFA3 and the quark masses. To 
change these, you must modify subroutine SSMSSM. By default, ALFA3=.12. 

All the mass spectrum and branching ratio calculations in ISASUSY are performed by a 
call to subroutine SSMSSM. Effective with version 7.23, the calling sequence is 

SUBROUTINE SSMSSM (XMG , XMU , XMHA , XTANB , XMQ 1 , XMDR , XMUR , 
$XML 1 , XMER , XMQ2 , XMSR , XMCR , XML2 , XMMR ,XMQ3, XMBR , XMTR , 
$XML3 , XMLR , XAT , XAB , XAL , XM1 , XM2 , XMT , I ALLOW , IMODEL) 

where the following are taken to be independent parameters: 
XMG = gluino mass 
XMU = n = SUSY Higgs mass 



XMHA 
XTANB 



—2 * mi of Baer et al. 
pseudo-scalar Higgs mass 
tan/3, ratio of vev's 
1/R (of old Baer- Tata notation). 
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XMQ1 = qi soft mass, 1st generation 

XMDR = d r mass, 1st generation 

XMUR = u r mass, 1st generation 

XML1 = £i mass, 1st generation 

XMER = e r mass, 1st generation 

XMQ2 = q~i soft mass, 2nd generation 

XMSR = s r mass, 2nd generation 

XMCR = c r mass, 2nd generation 

XML2 = li mass, 2nd generation 

XMMR = jj r mass, 2nd generation 

XMQ3 = qi soft mass, 3rd generation 

XMBR = b r mass, 3rd generation 

XMTR = t r mass, 3rd generation 

XML3 = £i mass, 3rd generation 

XMTR = f r mass, 3rd generation 

XAT = stop trilinear term A t 

XAB = sbottom trilinear term A^ 

XAL = stau trilinear term A T 

XM1 = U(l) gaugino mass 

= computed from XMG if I 1E19 

XM2 = SU(2) gaugino mass 

= computed from XMG if I 1E19 

XMT = top quark mass 
IALLOW = return flag 
IMODEL = 1 for SUGRA or MSSM 
= 2 for GMSB 
The variable IALLOW is returned: 
IALLOW = 1 if Z1SS is not LSP, otherwise 

All variables are of type REAL except IALLOW and IMODEL, which are INTEGER, and all 
masses are in GeV. The notation is taken to correspond to that of Haber and Kane, although 
the Tata Lagrangian is used internally. All other standard model parameters are hard wired 
in this subroutine; they are not obtained from the rest of ISAJET. The theoretically favored 
range of these parameters is 

100 < M(g) < 2000 GeV 
100 < M(q) < 2000 GeV 
100 < M(£) < 2000 GeV 
-1000 < /i < 1000 GeV 

1 < tan (3 < mt/rrib 
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M(t) « 175 GeV 
100 < M(A) < 2000 GeV 
M(tt),M(t r ) <M(q) 
M(b r ) ~ M(g) 
-1000 < A t < 1000 GeV 
-1000 <A b < 1000 GeV 

It is assumed that the lightest supersymmetric particle is the lightest neutralino Zi, the 
lighter stau fi, or the gravitino G in GMSB models. Some choices of the above parameters 
may violate this assumption, yielding a light chargino or light stop squark lighter than Z\. 
In such cases SSMSSM does not compute any branching ratios and returns IALLOW = I. 

SSMSSM does not check the parameters or resulting masses against existing experimental 
data. SSTEST provides a minimal test. This routine is called after SSMSSM by ISAJET 
and ISASUSY and prints suitable warning messages. 

SSMSSM first calculates the other SUSY masses and mixings and puts them in the 
common block /SSPAR/: 



c 


SUSY Darameters 






c 


AMGLSS 


_ 


gluino mass 


c 


AMULSS 


= 


un- left sauark mass 


c 


AMELSS 





lef t-selectron mass 


c 


AMERSS 




right-slepton mass 


c 


AMNiSS 




sneutrino mass for generation i 


c 


TW0M1 




Higgsino mass = - mu 


c 


RV2V1 




ratio v2/vl of vev's 


c 


AMTLSS , AMTRSS 




left, right stop masses 


c 


AMT1SS,AMT2SS 




light, heavy stop masses 


c 


AMBLSS , AMBRSS 




left, right sbottom masses 


c 


AMB1SS,AMB2SS 




light, heavy sbottom masses 


c 


AMLLSS , AMLRSS 




left, right stau masses 


c 


AML1SS,AML2SS 




light, heavy stau masses 


c 


AMZiSS 




signed mass of Zi 


c 


ZMIXSS 




Zi mixing matrix 


c 


AMWiSS 




signed Wi mass 


c 


GAMMAL , GAMMAR 




Wi left, right mixing angles 


c 


AMHL , AMHH , AMHA 




neutral Higgs hO, HO, AO masses 


c 


AMHC 




charged Higgs H+ mass 


c 


ALFAH 




Higgs mixing angle 


c 


AAT 




stop trilinear term 


c 


THETAT 




stop mixing angle 


c 


AAB 




sbottom trilinear term 


c 


THETAB 




sbottom mixing angle 


c 


AAL 




stau trilinear term 


c 


THETAL 




stau mixing angle 
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C AMGVSS = gravitino mass 

C MTQ = top mass at MSUSY 

C MBQ = bottom mass at MSUSY 

C MLQ = tau mass at MSUSY 

C FBMA = b- Yukawa at mA scale 

C VUQ = Hu vev at MSUSY 

C VDQ = Hd vev at MSUSY 



COMMON/SSPAR/AMGLSS , AMULSS , AMURSS , AMDLSS , AMDRSS , AMSLSS 
$ , AMSRSS , AMCLSS , AMCRSS , AMBLSS , AMBRSS , AMB1SS , AMB2SS 
$ , AMTLSS , AMTRSS , AMT1SS , AMT2SS , AMELSS , AMERSS , AMMLSS , AMMRSS 
$ , AMLLSS , AMLRSS , AML1SS , AML2SS , AMN1SS , AMN2SS , AMN3SS 
$ , TW0M1 , RV2V1 , AMZ1SS , AMZ2SS , AMZ3SS , AMZ4SS , ZMIXSS (4 ,4) 
$,AMW1SS,AMW2SS 

$ , GAMMAL , GAMMAR , AMHL , AMHH , AMHA , AMHC , ALFAH , AAT , THETAT 
$ , AAB , THETAB , AAL , THETAL , AMGVSS , MTQ , MBQ , MLQ , FBMA , 
$VUQ,VDQ 

REAL AMGLSS , AMULSS , AMURSS , AMDLSS , AMDRSS , AMSLSS 
$ , AMSRSS , AMCLSS , AMCRSS , AMBLSS , AMBRSS , AMB1SS , AMB2SS 
$ , AMTLSS , AMTRSS , AMT1SS , AMT2SS , AMELSS , AMERSS , AMMLSS , AMMRSS 
$ , AMLLSS , AMLRSS , AML1SS , AML2SS , AMN1SS , AMN2SS , AMN3SS 
$ , TW0M1 , RV2V1 , AMZ1SS , AMZ2SS , AMZ3SS , AMZ4SS , ZMIXSS 
$,AMW1SS,AMW2SS 

$ , GAMMAL , GAMMAR , AMHL , AMHH , AMHA , AMHC , ALFAH , AAT , THETAT 
$ , AAB , THETAB , AAL , THETAL , AMGVSS , MTQ , MBQ , MLQ , FBMA , VUQ , VDQ 
REAL AMZISS(4) 

EQUIVALENCE (AMZISS(l) ,AMZ1SS) 
SAVE /SSPAR/ 

It then calculates the widths and branching ratios and puts them in the common block 
/SSMODE/: 



C MXSS = maximum number of modes 

C NSSMOD = number of modes 

C ISSMOD = initial particle 

C JSSMOD = final particles 

C GSSMOD = width 

C BSSMOD = branching ratio 

C MSSMOD = decay matrix element pointer 

C LSSMOD = logical flag used internally by SSME3 



INTEGER MXSS 
PARAMETER (MXSS=1000) 

COMMON/SSMODE/NSSMOD , ISSMOD (MXSS) , JSSMOD (5 , MXSS) , GSSMOD (MXSS) 
$, BSSMOD (MXSS) , MSSMOD (MXSS) , LSSMOD 
INTEGER NSSMOD , ISSMOD , JSSMOD , MSSMOD 
REAL GSSMOD, BSSMOD 
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LOGICAL LSSMOD 
SAVE /SSMODE/ 



Decay modes for a given particle are not necessarily adjacent in this common block. Note 
that the branching ratio calculations use the full matrix elements, which in general will give 
nonuniform distributions in phase space, but this information is not saved in /SSMODE/. 
In particular, the decays H-^Z + Z*^Z + f + f give no indication that the // mass is 
strongly peaked near the upper limit. 

All IDENT codes are defined by parameter statements in the PATCHY keep sequence 
SSTYPE: 

C SM ident code definitions. These are standard ISAJET but 

C can be changed. 

INTEGER IDUP , IDDN , IDST , IDCH , IDBT , IDTP 

INTEGER IDNE , IDE , IDNM , IDMU , IDNT , IDTAU 

INTEGER IDGL , IDGM , IDW , IDZ , IDH 

PARAMETER (IDUP=1 , IDDN=2 , IDST=3 , IDCH=4 , IDBT=5 , IDTP=6) 

PARAMETER (IDNE=11 , IDE=12 , IDNM=13 , IDMU=14 , IDNT=15 , IDTAU=16) 

PARAMETER (IDGL=9 , IDGM=10 , IDW=80 , IDZ=90 , IDH=81) 
C SUSY ident code definitions. They are chosen to be similar 

C to those in versions < 6.50 but may be changed. 

INTEGER ISUPL , ISDNL , ISSTL , ISCHL , ISBT1 , ISTP1 

INTEGER ISNEL , ISEL , ISNML , ISMUL , ISNTL , ISTAU1 

INTEGER ISUPR, ISDNR, ISSTR, ISCHR, ISBT2, ISTP2 

INTEGER ISNER, ISER, ISNMR, ISMUR, ISNTR, ISTAU2 

INTEGER ISZ1 , ISZ2 , ISZ3 , ISZ4 , ISW1 , ISW2 , ISGL 

INTEGER ISHL,ISHH,ISHA,ISHC 

INTEGER ISGRAV 

INTEGER IDTAUL , IDTAUR 

PARAMETER (ISUPL=21 , ISDNL=22 , ISSTL=23 , ISCHL=24 , ISBT1=25 , ISTP1=26) 
PARAMETER (ISNEL=31 , ISEL=32 , ISNML=33 , ISMUL=34 , ISNTL=35 , ISTAU1=36) 
PARAMETER (ISUPR=41 , ISDNR=42 , ISSTR=43 , ISCHR=44 , ISBT2=45 , ISTP2=46) 
PARAMETER (ISNER=51 , ISER=52 , ISNMR=53 , ISMUR=54 , ISNTR=55 , ISTAU2=56) 
PARAMETER (ISGL=29) 

PARAMETER (ISZ1=30 , ISZ2=40 , ISZ3=50 , ISZ4=60 , ISW1=39 , ISW2=49) 

PARAMETER (ISHL=82 , ISHH=83 , ISHA=84, ISHC=86) 

PARAMETER (ISGRAV=91) 

PARAMETER (IDTAUL=10016 , IDTAUR=20016) 

These are based on standard ISAJET but can be changed to interface with other generators. 
Since masses except the t mass are hard wired, one should check the kinematics for any 
decay before using it with possibly different masses. 

Instead of specifying all the SUSY parameters at the electroweak scale using the MSSMi 
commands, one can instead use the SUGRA parameter to specify in the minimal supergravity 
framework the common scalar mass m , the common gaugino mass m^, and the soft trilinear 
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SUSY breaking parameter A at the GUT scale, the ratio tan/5 of Higgs vacuum expectation 
values at the electroweak scale, and sgn /i, the sign of the Higgsino mass term. The NUSUGi 
keywords allow one to break the assumption of universality in various ways. NUSUGI sets the 
gaugino masses; NUSUG2 sets the A terms; NUSUG3 sets the Higgs masses; NUSUG4 sets the first 
generation squark and slepton masses; and NUSUG5 sets the third generation masses. The 
keyword SSBCSC can be used to specify an alternative scale (i.e., not the coupling constant 
unification scale) for the RGE boundary conditions. 

The renormalization group equations now include all the two-loop terms for both gauge 
and Yukawa couplings and the possible contributions from right-handed neutrinos. These 
equations are solved iteratively using Runge-Kutta numerical integration to determine the 
weak scale parameters from the GUT scale ones: 

1. The 2-loop RGE's for the gauge, Yukawa, and soft breaking terms are run from the 
weak scale Mz up to the GUT scale (where oc± — 0,2) or other high scale defined by 
the model, taking all thresholds into account. 

2. The GUT scale boundary conditions are imposed, and the RGE's are run back to Mz, 
again taking thresholds into account. 

3. The masses of the SUSY particles and the values of the soft breaking parameters B 
and mu needed for radiative symmetry are computed, e.g. 



4. The 1-loop radiative corrections are computed. 

5. The process is then iterated until stable results are obtained. 

This is essentially identical to the procedure used by several other groups. Other possible 
constraints such as b-r unification and limits on proton decay have not been included. 

An alternative to the SUGRA model is the Gauge Mediated SUSY Breaking (GMSB) 
model of Dine and Nelson, Phys. Rev. D48, 1277 (1973); Dine, Nelson, Nir, and Shirman, 
Phys. Rev. D53, 2658 (1996). In this model SUSY is broken dynamically and communicated 
to the MSSM through messenger fields at a messenger mass scale M m much less than the 
Planck scale. If the messenger fields are in complete representations of SU(5), then the 
unification of couplings suggested by the LEP data is preserved. The simplest model has 
a single 5 + 5 messenger sector with a mass M m and and a SUSY-breaking VEV F m of its 
auxiliary field F. Gauginos get masses from one-loop graphs proportional to A m = F m /M m 
times the appropriate gauge coupling af, sfermions get squared- masses from two- loop graphs 
proportional to A m times the square of the appropriate ctj. If there are N 5 messenger fields, 
the gaugino masses and sfermion masses- squared each contain a factor of N 5 . 

The parameters of the GMSB model implemented in ISAJET are 



Ml 



L - M 2 H2 tan 2 f3 
tan 2 - 1 



Ml/2 




• A 



m 



F m /M m : the scale of SUSY breaking, typically 10-100 TeV; 
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• M m > A m : the messenger mass scale, at which the boundary conditions for the renor- 
malization group equations are imposed; 

• N 5 : the equivalent number of 5 + 5 messenger fields. 

• tan/5: the ratio of Higgs vacuum expectation values at the electroweak scale; 

• sgn/i = ±1: the sign of the Higgsino mass term; 

• Cgrav > 1: the ratio of the gravitino mass to the value it would have had if the only 
SUSY breaking scale were F m . 

The solution of the renormalization group equations is essentially the same as for SUGRA; 
only the boundary conditions are changed. In particular it is assumed that electroweak 
symmetry is broken radiatively by the top Yukawa coupling. 

In GMSB models the lightest SUSY particle is always the nearly massless gravitino G. 
The phenomenology depends on the nature of the next lightest SUSY particle (NLSP) and 
on its lifetime to decay to a gravitino. The NLSP can be either a neutralino Xi or a slepton 
f\. Its lifetime depends on the gravitino mass, which is determined by the scale of SUSY 
breaking not just in the messenger sector but also in any other hidden sector. If this is set 
by the messenger scale F m , i.e., if C grav ~ 1, then this lifetime is generally short. However, 
if the messenger SUSY breaking scale F rn is related by a small coupling constant to a much 
larger SUSY breaking scale F b , then C grav ^> 1 and the NLSP can be long-lived. The correct 
scale is not known, so C grav should be treated as an arbitrary parameter. More complicated 
GMSB models may be run by using the GMSB2 keyword. 

Patch ISASSRUN of ISAJET provides a main program SSRUN and some utility programs 
to produce human readable output. These utilities must be rewritten if the IDENT codes in 
/SSTYPE/ are modified. To create the stand-alone version of ISASUSY with SSRUN, run 
YPATCHY on isajet.car with the following cradle (with & replaced by +): 

&USE,* ISASUSY. Select all code 

&USE,NOCERN. No CERN Library 

&USE,IMPNONE. Use IMPLICIT NONE 

&EXE. Write everything to ASM 

&PAM,T=C. Read PAM file 

&QUIT. Quit 

Compile, link, and run the resulting program, and follow the prompts for input. Patch 
ISASSRUN also contains a main program SUGRUN that reads the minimal SUGRA, non- 
universal SUGRA, or GMSB parameters, solves the renormalization group equations, and 
calculates the masses and branching ratios. To create the stand-alone version of ISASUGRA, 
run YPATCHY with the following cradle: 

&USE, * ISASUGRA . Select all code 

&USE,NOCERN. No CERN Library 

&USE,IMPNONE. Use IMPLICIT NONE 

&EXE. Write everything to ASM 
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&PAM. 
&QUIT. 



Read PAM file 
Quit 



The documentation for ISASUSY and ISASUGRA is included with that for ISAJET. 

ISASUSY is written in ANSI standard Fortran 77 except that IMPLICIT NONE is used 
if +USE,IMPNONE is selected in the Patchy cradle. All variables are explicitly typed, and 
variables starting with I,J,K,L,M,N are not necessarily integers. All external names such as 
the names of subroutines and common blocks start with the letters SS. Most calculations 
are done in double precision. If +USE,NOCERN is selected in the Patchy cradle, then the 
Cernlib routines EISRS1 and its auxiliaries to calculate the eigenvalues of a real symmetric 
matrix and DDILOG to calculate the dilogarithm function are included. Hence it is not 
necessary to link with Cernlib. 

The physics assumptions and details of incorporating the Minimal Supersymmetric Model 
into ISAJET have appeared in a conference proceedings entitled "Simulating Supersymmetry 
with ISAJET 7.0/ISASUSY 1.0" by H. Baer, F. Paige, S. Protopopescu and X. Tata; this 
has appeared in the proceedings of the workshop on Physics at Current Accelerators and 
Supercolliders, ed. J. Hewett, A. White and D. Zeppenfeld, (Argonne National Laboratory, 
1993). Detailed references may be found therein. Users wishing to cite an appropriate source 
should, however, cite the most recent ISAJET manual, e.g. hep-ph/0001086 (1999). 
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13 Changes in Recent Versions 



This section contains a record of changes in recently released versions of ISAJET, taken 
from the memoranda distributed to users. Note that the released version numbers are not 
necessarily consecutive. 

13.1 Version 7.69, August 2003 

The complete set of 1-loop radiative corrections to sparticle masses has been included in the 
ISAJET SUSY spectrum calculation, according to the formulae given by Pierce et al, Nucl. 
Phys. B491, 3 (1997). Previous versions of ISAJET had included the logarithmic 1-loop 
corrections by freezing out soft mass parameters at a scale equal to their mass. The current 
calculation includes all finite corrections as well. Crucial contributions to the encoding of 
these expressions were made by Tadas Krupovnickas. We have also adjusted slightly the input 
DR gauge couplings at Mz to coincide with measured central values. Logarithmic threshold 
corrections to gauge and Yukawa couplings are accounted for as in previous ISAJET versions 
by changing the beta functions in the RGEs as various soft mass thresholds are passed. We 
have also enlarged the set of SUSY RGEs to include running of the Higgs field vevs. The 
net effect of these changes is to modify various sparticle masses by typically 1-5% from the 
predictions of ISAJET 7.64. 

The size of the small mass splitting between xf and Xi i s important for AMSB phe- 
nomenology. Rather than extracting this splitting from the general result, for the AMSB 
model we calculate it from the tree and the 1-loop vector boson loop graphs [taken from 
Cheng, Dobrescu and Matchev, Nucl. Phys. B543, 47, (1999)] and use it to set the xf mass. 

We have modified the RGE solution to use the previous SUSY masses to compute fi 
before computing the new SUSY masses and from them the loop corrections to fx. This 
seems to improve the convergence, e.g., in the focus point, or hyperbolic branch region, of 
the mSUGRA model. The RGE solution still does not converge properly in isolated regions. 

We have also corrected a sign mistake in the A-term boundary conditions for the anomaly- 
mediated SUSY breaking model. 

13.2 Version 7.64, September 2002 

The talk by Sabine Kraml at SUSY02 stressed the sensitivity of the allowed region of radiative 
electroweak symmetry breaking to fine details of the calculation. We have reexamined the 
issue and corrected several problems. A coding error in the Passarino-Veltman function B\, 
SSB1 has been fixed. The requirement fx 2 > is not imposed until after the solution of the 
renormalization group equations has stabilized. The value of m^imz) has been updated to 
2.83 GeV. A running gluino mass has been used in the top and bottom self-energy, along 
with a s evaluated at a higher scale. Since the parameters of the Higgs potential vary rapidly 
near the weak scale, the convergence requirements on these have been loosened, while the 
requirements on the other parameters have been somewhat tightened. The combined effect 
of these changes is to shift the boundary for radiative electroweak symmetry breaking to 
higher scalar masses. 
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Stephan Lammel found that the matrix element used to generate g — ► xtlQ was missing 
some poles, although the branching ratio was correct. This has been fixed. 
The radiative decays Xi ~* X°7 have been included. 

13.3 Version 7.63, April 2002 

The SUSY mass calculations have been improved, especially for M A in terms of other SUSY 
parameters, by using the MSSM Yukawa couplings from the renormalization group equations. 
The numerical precision of the solution to the SUSY renormalization group equations has 
also been improved; this should give better stability near the boundaries of the allowed 
regions. The complete 1-loop self-energies for the t, b, and r have been included from Pierce, 
Bagger, Matchev, and Zhang, Nucl. Phys. B491, 3 (1997). Finally, a number of bugs have 
been fixed, including one in the r decay of t quarks. 

13.4 Version 7.58, August 2001 

The CTEQ5L parton distributions have been added and made the default. 

Keywords NOB and NOTAU have been added to turn off B and r decays so that an external 
decay package such as QQ or TAUOLA can be used. To preserve the polarization information 
for r's, separate (non-standard) IDENT codes for tl and tr are used for NOTAU. The user 
must provide an appropriate interface. 

The RANLUX random number generator has been added as a compile time option. It has 
a very long period, and any 32-bit integer seed gives an independent sequence. If RANLUX is 
used, the keyword SEED takes the integer seed plus two additional integers that are normally 
zero but can be used to restart the generator. See the CERN Program Library writeups for 
more information. 

Right-handed neutrinos are now included if the keyword SUGRHN is used. The user must 
input the 3rd generation neutrino mass (at scale M z ), the intermediate scale right handed 
neutrino Majorana mass M N , and the soft SUSY-breaking masses A n and at the GUT 
scale. Then the neutrino Yukawa coupling is computed in the simple see-saw model, and 
renormalization group evolution includes these effects between Mqut and Mjv- 

The decays g — > W(ud have been updated to include non-degenerate squark masses. The 
arbitrary width for t — > Z\C used previously has been replaced by the calculated value. 
DR masses are used consistently. Yukawa couplings in the the SUGRA routine are now 
calculated in the DR regularization scheme to be consistent with two loop renormalization 
group evolution. 

In solving the SUSY renormalization group equations, the requirement of good elec- 
troweak symmetry breaking is imposed only at the end. Previously a point could be rejected 
if there was no symmetry breaking even in the initial iteration with a truncated set of 
equations. 

The sign of the A term in the AMSB model has been corrected. 
The Standard Model process e + e~ — > ZH was missing and has been added. 
Function ITRANS has been updated to reflect the current PDG particle codes for SUSY 
particles. 



88 



13.5 Version 7.51, May 2000 

Several improvements in the SUSY RGE's have been made. All two-loop terms including 
both gauge and Yukawa couplings and the contributions from right-handed neutrinos are 
now included. There is a new keyword SSBCSC to specify a scale other than the GUT scale 
for the RGE boundary conditions. 

The process Z + 7 is now included in WPAIR. (This was omitted because it has no contri- 
bution from triple gauge boson couplings.) 

An incorrect type declaration produced unphysical results for beamsstrahlung on some 
computers. This has been fixed. While the bug is serious for e + e~ with the EEBEAM option, 
it has no effect on other processes. Some other minor bugs have also been fixed. 

13.6 Version 7.47, December 1999 

There are several improvements in the treatment of supersymmetry. The Anomaly Mediated 
SUSY Breaking model of of Randall and Sundrum and of Gherghetta, Giudice, and Wells 
(hep-ph/9904378) has been added. The parameters of the model are a universal scalar mass 
m at the GUT scale, a gravitino mass m 3 / 2 , and the usual tan/3 and sgn/i. These are set 
by the AMSB keyword. The renormalization group equations have been extended to include 
two-loop Yukawa terms and right-handed sneutrinos (with default masses above the Planck 
scale). The u R play a role in the evolution for the inverted hierarchy models of Bagger, 
Feng, and Polonsky, hep-ph/9905292. SUSY loop corrections to Yukawa couplings have 
been incorporated in the SUSY mass calculations. 

The Helas library of Murayama, Watanabe, and Hagiwara has been incorporated together 
with a simple multi-body phase space generator. This makes it possible to use code generated 
by MadGraph to produce multi-body hard scattering processes. As a first example, a ZJJ 
process that generates Z + 2 jets has been added, with the Z treated as a narrow resonance. 
Additional processes may be added in future releases. 

A new EXTRADIM process has been added to generate Kaluza-Klein graviton production 
in association with a jet or photon in models with extra dimensions at the TeV scale. The 
cross sections are from G.F. Giudice et al., hep-ph/9811291. We thank I. Hinchliffe and L. 
Vacavant for providing this. 

A number of bugs have been fixed, including in particular one in the decay Wi — > Z^rv. 

13.7 Version 7.44, April 1999 

A serious bug introduced in Version 7.42 that could lead to matrix elements being stored for 
the wrong mode has been corrected. Some sign errors in the matrix elements for gaugino 
decays have also been corrected. 

13.8 Version 7.42, January 1999 

Beginning with this version, matrix elements are taken into account in the event generator 
as well as in the calculation of decay widths for MSSM three-body decays of the form 
A — > Bff, where A and B are gluinos, charginos, or neutralinos. This is implemented 
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by having ISASUSY save the poles and their couplings when calculating the decay width 
and then using these to reconstruct the matrix element. Other three-body decays may be 
included in the future. Decays selected with FORCE use the appropriate matrix elements. 

As part of the changes to implement these matrix elements, the format of the decay table 
has changed. It now starts with a header line; if this does not match the internal version, 
then a warning is printed. The decay table now includes an index MELEM that specifies 
the matrix element to be used for all processes. This is also used for FORCE decays and is 
printed on the run listing for them. SUSY 3-body decays have internally generated negative 
values of MELEM. 

This version also includes both initial state radiation and beamstrahlung for e + e~ in- 
teractions. For initial state radiation (bremsstrahlung) , if the EEBREM keyword is selected, 
an electron structure function will be used. For a convolution of both bremsstrahlung and 
beamstrahlung, the keyword EEBEAM must be used, with appropriate inputs (see documen- 
tation). 

13.9 Version 7.40, October 1998 

A new process WHIGGS generates W ± + H and Z + H events for both the Standard Model 
and SUSY models and also Higgs pair production for SUSY models. The types and W 
decay modes are selected with JETTYPE and WMODE as for WPAIR events. This process 
is of particular interest for producing fairly light Higgs bosons at the Tevatron. See the 
documentation for more details. 

Some non-minimal GMSB models can be generated using a new keyword GMSB2. The 
optional parameters are an extra factor between the gaugino and scalar masses, shifts in the 
Higgs masses, a D-term proportional to hypercharge, and independent numbers of messenger 
fields for the three gauge groups. The documentation gives more details and references. 

The default for SUGRA models has been changed to use a s (M z ) = 0.118, the experimen- 
tal value. This means that the couplings do not exactly unify at the GUT scale, presumably 
because of the effects of heavy particles. The keyword AL3UNI can be used to select exact 
unification, which produces too large a value for a s (Mz)- 

A number of three-body slepton decays that occur through left-right mixing are now 
included. These are obviously small but might compete with gravitino decays. In particular, 
a decay like jiR — > f\vv might lead to a wrong momentum measurement in the muon system. 
So far we have found no case in which this is probable. 

The new release also includes a separate Unix tar file mcpp.tar containing C++ code to 
read a standard ISAJET output file and copy all the information into C++ classes. The tar 
file contains makefiles for Software Release Tools, documentation, and examples as well as 
the code. 

13.10 Version 7.37, April 1998 

Version 7.37 incorporates Gauge Mediated SUSY Breaking models for the first time. In 
these models, SUSY is broken in a hidden sector at a relatively low scale, and the masses 
of the MSSM fields are then produced through ordinary gauge interactions with messenger 
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fields. The parameters of the GMSB model in ISAJET are M m , the messenger mass scale; 
A m = F m /M m , where F m is the SUSY breaking scale in the messenger sector; N 5 , the 
number of messenger fields; the usual tan/5 and sgn/x; and C grav > 1, a factor which scales 
the gravitino mass and hence the lifetime for the lightest MSSM particle to decay into it. 

GMSB models have a light gravitino G as the lightest SUSY particle. The phenomenology 
of the model depends mainly on the nature of the next lightest SUSY particle, a x? or a 
which changes with the number N 5 of messengers. The phenomenology also depends on the 
lifetime for the Xi ~^ ^7 or ^1 — * Gt decay; this lifetime can be short or very long. All the 
relevant decays are included except for \x — > vvf\, which is very suppressed. 

The keyword MGVTNO allows the user to independently input a gravitino gravitino 
mass for the MSSM option. This allows studies of SUGRA (or other types) of models where 
the gravitino is the LSP. 

Version 7.37 also contains an extension of the SUGRA model with a variety of non- 
universal gaugino and sfermion masses and A terms at the GUT scale. This makes it possible 
to study, for example, how well the SUGRA assumptions can be tested. 

Two significant bugs have also been corrected. The decay modes for B* mesons were 
missing from the decay table since Version 7.29 and have been restored. A sign error in the 
interference term for chargino production has been corrected, leading to a larger chargino 
pair cross section at the Tevatron. 

13.11 Version 7.32, November 1997 

This version makes several corrections in various chargino and neutralino widths, thus chang- 
ing the branching ratios for large tan/3. For \% f° r example, the Xi^b branching ratio is 
decreased significantly, while the Xi T+T ~ one is increased. Thus the SUGRA phenomenology 
for tan (3 ~ 30 is modified substantially. 

The new version also fixes a few bugs, including a possible numerical precision problem 
in the Drell-Yan process at high mass and qx- It also includes a missing routine for the Zebra 
interface. 

13.12 Version 7.31, August 1997 

Version fixes a couple of bugs in Version 7.29. In particular, the JETTYPE selection did not 
work correctly for supersymmetric Higgs bosons, and there was an error in the interactive 
interface for MSSM input. Since these could lead to incorrect results, users should replace 
the old version. We thank Art Kreymer for finding these problems. 

Since top quarks decay before they have time to hadronize, they are now put directly 
onto the particle list. Top hadrons (tu, td, etc.) no longer appear, and FORCE should be 
used directly for the top quark, i.e. 

FORCE 

6,11,-12,5/ 

The documentation has been converted to LaTeX. Run either LaTeX 2.09 or LaTeX 2e 
three times to resolve all the forward references. Either US (8.5x11 inch) or A4 size paper 
can be used. 
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13.13 Version 7.30, July 1997 



This version fixes a couple of bugs in the previous version. In particular, the JETTYPE 
selection did not work correctly for supersymmetric Higgs bosons, and there was an error in 
the interactive interface for MSSM input. Since these could lead to incorrect results, users 
should replace the old version. We thank Art Kreymer for finding these problems. 

Since top quarks decay before they have time to hadronize, they are now put directly 
onto the particle list. Top hadrons (tu, tud, etc.) no longer appear, and FORCE should be 
used directly for the top quark, i.e. 

FORCE 

6,11,-12,5/ 

The documentation has been converted to M^X. Run either DTj^X 2.09 or I^TgX 2e 
three times to resolve all the forward references. Either US (8.5 x 11 inch) or A4 size paper 
can be used. 

13.14 Version 7.29, May 1997 

While the previous version was applicable for large as well as small tan/3, it did contain 
approximations for the 3-body decays g — > tbWi, Zi — > hhZ r ; : and Wi — > ruZj. The 
complete tree-level calculations for three body decays of the gluino, chargino and neutralino, 
with all Yukawa couplings and mixings, have now been included (thanks mainly to M. Drees). 
We have compared our branching ratios with those calculated by A. Bartl and collaborators; 
the agreement is generally good. 

The decay patterns of gluinos, charginos and neutralinos may differ from previous ex- 
pectations if tan (3 is large. In particular, decays into r's and fo's are often enhanced, while 
decays into e's and /z's are reduced. It could be important for experiments to study new 
types of signatures, since the cross sections for conventional signatures may be considerably 
reduced. 

We have also corrected several bugs, including a fairly serious one in the selection of jet 
types for SUSY Higgs. We thank A. Kreymer for pointing this out to us. 

13.15 Version 7.27, January 1997 

The new version contains substantial improvements in the treatment of the Minimal Su- 
persymmetric Standard Model (MSSM) and the SUGRA model. The squarks of the first 
two generations are no longer assumed to be degenerate. The mass splittings and all the 
two-body decay modes are now correctly calculated for large tan/5. While there are still 
some approximations for three-body modes, ISAJET is now usable for the whole range 
1 ^ tan (3 ^ M t /M b . The most interesting new feature for large tan/5 is that third genera- 
tion modes can be strongly enhanced or even completely dominant. 

To accomodate these changes it was necessary to change the MSSM input parameters. To 
avoid confusion, the MSSM keywords have been renamed MSSM[A-C] instead of MSSM[l-3], 
and the order of the parameters has been changed. See the input section of the manual for 
details. 
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Treatment of the MSSM Higgs sector has also been improved. In the renormalization 
group equations the Higgs couplings are frozen at a higher scale, Q = \J M(t L )M{t R ). Run- 
ning t, b and r masses evaluated at that scale are used to reproduce the dominant 2-loop 
effects. There is some sensitivity to the choice of Q; our choice seems to give fairly stable 
results over a wide range of parameters and reasonable agreement with other calculations. 
In particular, the resulting light Higgs masses are significantly lower than those from Ver- 
sion 7.22. 

The default parton distributions have been updated to CTEQ3L. A bug in the PDFLIB 
interface and other minor bugs have been fixed. 

13.16 Version 7.22, July 1996 

The new version fixes errors in b — > Wt and in some t decays and Higgs decays. It also 
contains a new decay table with updated r, c, and b decays, based loosely on the QQ 
decay package from CLEO. The updated decays are less detailed than the full CLEO QQ 
program but an improvement over what existed before. The new decays involve a number of 
additional resonances, including / (980), ai(1260), / 2 (1270), 7^(1270), 7^(1400), X 2 *(1430), 
Xci,2,3, and ip{2S), so users may have to change their interface routines. 
A number of other small bugs have been corrected. 

13.17 Version 7.20, June 1996 

The new version corrects both errors introduced in Version 7.19 and longstanding errors in 
the final state QCD shower algorithm. It also includes the top mass in the cross sections for 
gb — > Wt and gt — > Zt. When the t mass is taken into account, the process gt — > Wb can 
have a pole in the physical region, so it has been removed; see the documentation for more 
discussion. 

Steve Tether recently pointed out to us that the anomalous dimension for the q — > qg 
branching used in the final state QCD branching algorithm was incorrect. In investigating 
this we found an additional error, a missing factor of 1/3 in the g — > qq branching. The 
first error produces a small but non-negligible underestimate of gluon radiation from quarks. 
The second overestimates quark pair production from gluons by about a factor of 3. In 
particular, this means that backgrounds from heavy quarks Q coming from g — > QQ have 
been overestimated. 

The new version also allows the user to set arbitrary masses for the U(l) and SU{2) 
gaugino mases in the MSSM rather than deriving these from the gluino mass using grand 
unification. This could be useful in studying one of the SUSY interpretations of a CDF 
ee77^T event recently suggested by Ambrosanio, Kane, Kribs, Martin and Mrenna. Note, 
however, that radiative decay are not included, although the user can force them and multiply 
by the appropriate branching ratios calculated by Haber and Wyler, Nucl. Phys. B323, 267 
(1989). No explicit provision for the decay Z\ — > G7 of the lightest zino into a gravitino or 
goldstino and a photon has been made, but forcing the decay Z\ — > 1/7 has the same effect 
for any collider detector. 

A number of other minor bugs have also been corrected. 
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13.18 Version 7.16, October 1995 



The new version includes e + e~ cross sections for both SUSY and Standard Model particles 
with polarized beams. The e" and e + polarizations are specified with a new keyword EPOL. 
Polarization appears to be quite useful in studying SUSY particles at an e + e~ collider. 

The new release also includes some bug fixes for pp reactions, so you should upgrade even 
if you do not plan to use the polarized e + e~ cross sections. 

13.19 Version 7.13, September 1994 

Version 7.13 of ISAJET fixes a bug that we introduced in the recently released 7.11 and 
another bug in g — > qq. We felt it was essential to fix these bugs despite the proliferation of 
versions. 

The new version includes the cross sections for the e + e~ production of squarks, sleptons, 
gauginos, and Higgs bosons in Minimal Supersymmetric Standard Model (MSSM) or the 
minimal supergravity (SUGRA) model, including the effects of cascade decays. To generate 
such events, select the E+E- reaction type and either SUGRA or MSSM, e.g., 

SAMPLE E+E- JOB 
300. ,50000,10,100/ 
E+E- 
SUGRA 

100,100,0,2,-1/ 

TMASS 

170,-1,1/ 

END 

STOP 

The effects of spin correlations in the production and decay, e.g., in e + e~ — > W'+Wf, are 
not included. 

It should be noted that the Standard Model e + e~ generator in ISAJET does not include 
Bhabba scattering or W + W~ and Z°Z° production. Also, its hadronization model is cruder 
than that available in some other generators. 

13.20 Version 7.11, September 1994 

The new version includes the cross sections for the e + e~ production of squarks, sleptons, 
gauginos, and Higgs bosons in Minimal Supersymmetric Standard Model (MSSM) or the 
minimal supergravity (SUGRA) model including the effects of cascade decays. To generate 
such events, select the E+E- reaction type and either SUGRA or MSSM, e.g., 

SAMPLE E+E- JOB 
300. ,50000,10,100/ 
E+E- 
SUGRA 

100,100,0,2,-1/ 
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TMASS 

170,-1,1/ 

END 

STOP 

The effects of spin correlations in the production and decay, e.g., in e + e~ — * Wj^Wf, are 
not included. 

It should be noted that the Standard Model e + e~ generator in ISAJET does not include 
Bhabba scattering or W + W~ and Z°Z° production. Also, its hadronization model is cruder 
than that available in some other generators. 

13.21 Version 7.10, July 1994 

This version adds a new option that solves the renormalization group equations to calculate 
the Minimal Supersymmetric Standard Model (MSSM) parameters in the minimal super- 
gravity (SUGRA) model, assuming only that the low energy theory has the minimal particle 
content, that electroweak symmetry is radiatively broken, and that R-parity is conserved. 
The minimal SUGRA model has just four parameters, which are taken to be the common 
scalar mass too, the common gaugino mass m^, the common trilinear SUSY breaking term 
Aq, all defined at the GUT scale, and tan /3; the sign of [i must also be given. The renormal- 
ization group equations are solved iteratively using Runge-Kutta integration including the 
correct thresholds. This program can be used either alone or as part of the event generator. 
In the latter case, the parameters are specified using 

SUGRA 

mo, toi/2, A , tan/3, sgn/i 

While the SUGRA option is less general than the MSSM, it is theoretically attractive and 
provides a much more managable parameter space. 

In addition there have been a number of improvements and bug fixes. An occasional 
infinite loop in the minimum bias generator has been fixed. A few SUSY cross sections 
and decay modes and the JETTYPE flags for SUSY particles have been corrected. The 
treatment of B baryons has been improved somewhat. 
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